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The  main  body  of  this  report  is  being  presented  by  Mr.  Jerone  J. 
Vorachek  (Project  Engineer)  and  Mr.  George  R,  Doyle,  Jr.  The 
equations  of  motion  were  derived  by  Mr.  Doyle  who  would  uve  to 
thank  Mr.  Bernard  Burzlaff  and  Mr.  Nicholas  Odrey  for  offering 
suggestions  and  reviewing  the  mathematical  approach  and  derivation 
of  the  equations.  This  derivation  is  presented  in  Appendix  A. 
Appendix  B  contains  an  explanation  of  the  computer  program  used  to 
obtain  the  results.  This  program  was  written  by  Mr.  Doyle. 
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The  desire  to  place  a  payload  at  altitudes  of  100,000  ft  and  to  be  able 
to  retrieve  it  prompted  a  study  of  the  ascent  and  descent  dynamics  of  a 
balloon-tether  system.  The  basic  incentive  of  this  study  is  to  establish 
the  feasibility  of  deploying  and  retrieving  high  altitude  tethered  balloons 
through  high  wind  velocity  regions  at  intermediate  altitudes.  It  was  necessary 
to  establish  the  technical  ability  to  analyze  such  a  system,  and  to  develop  a 
computer  program  which  would  define  the  dynamics  of  the  tether  and  balloon  by 
integrating  the  equations  of  motion.  Lagrange's  equation  was  used  to  derive 
the  equations  of  motion;  Runge-Kutta  numerical  integration  was  used  to  solve 
the  equations  of  motion;  and  the  computer  program  was  written  in  Fortran  IV 
for  the  IBM  360, 

The  system  is  assumed  to  be  comprised  of  rigid  bodies  connected  by  frictionless 
hinges.  The  tether  is  simulated  by  "N"  straight  links.  All  of  the  links 
are  the  same  length  and  grow  at  the  same  rate,  but  do  not  have  the  same  weight 
or  aerodynamic  reference  area.  For  this  preliminary  study,  the  balloon 
is  considered  to  be  a  thin  spherical  shell  which  expands  as  it  rises  in 
proportion  to  density  ratio  change  with  altitude. 

The  results  of  a  limited  number  of  computer  runs  is  being  presented  not  as  an 
extensive  study  of  balloon-tether  systems,  but  rather  as  a  demonstration 
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8”Ch  S7BteM  fea»lMe  —  •— »  conditions.  Also  present*  1.  an 
understand  of  the  Abilities  and  Station.  of  *.  softer  progRM 

“  the  s.stan.  The  Mjor  enphasi.  is  piaced  on  ascent  and  desoent 

^  a  Sumner  I  wind  prod.  A  more  Tinted  effort  is  aiso  Resented  on 
ascent  and  descent  of  a  balloon  in  a  Winter  I  wind  proftle. 

ngures  1  and  2  show  the  Su-er  I  and  Winter  I  wind  profiles  as  the  counter 

interprets  them.  Sumner  I  75  percent  and  Winter  I  75  n«ro  + 

,  ,  ter  I  75  percent  wind  profiles  may 

be  found  in  Reference  3.  Summer  I  90  nercant  oc 

percent,  95  percent  and  Winter  I  25  Percent 
come  from  Reference  h.  Percent 
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SECTION  II  -  RESORTS  OF  ASCENT  AND  DESCENT  FEASIBIHTY  STUDY 


The  feasibility  study  was  conducted  using  two  balloon-tether  system  - 

one  suitable  for  Sumer  I  wind  conditions  and  one  for  Winter  I  wind  conditions, 

A  description  of  the  balloon-tether  systems  is  given  in  Appendix  B-2  and  B-6. 

The  design  of  the  systems  is  further  discussed  in  Reference  5,  The  basic  criteria 
of  each  system  is  that  the  balloon  support  a  cable  (factor  of  safety  equal  2) 
at  100,000  ft  in  a  75  percent  wind  profile.  This  defines  static  conditions. 

It  is  Uie  purpose  of  this  study  to  determine  whether  or  not  the  same  system  can 
satisfy  dynamic  conditions  of  ascent  and  descent.  In  general,  it  has  been 
found  that  soiob  conditions  can  be  satisfied  and  some  can  not,  but  this  study 
does  not  attempt  to  define  limiting  conditions.  Effects  of  wind  profile, 
winching  rate  and  balloon  size  are  established. 

In  order  to  show  that  a  zero  drag  tether  does  not  adequately  simulate  the 
balloon  trajectory.  Figure  3  is  presented.  An  obvious  divergence  can  be  seen 
above  20,000  ft.  At  h0,000  ft,  the  divergence  becomes  more  pronounced.  This 
can  be  attributed  to  the  fact  that  the  region  of  high  dynamic  pressure  is 
between  35,000  and  1*0,000  ft.  In  one  trajectory,  the  system  still  experiences 
the  high  dynamic  pressure  through  the  tether  even  though  the  balloon  is  approaching 
a  lower  dynamic  pressure  region.  In  the  other,  the  system  experiences  only  the 
dynamic  pressure  at  the  altitude  of  the  balloon  causing  it  to  start  to  drift  back 
to  the  vertical. 
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Figures  U,  5,  6,  and  7  demonstrate  the  effects  of  pay-out  rate.  It  is  realized 
that  a  rate  of  3000  ft/min  is  pushing  the  w  state-of-the-art”  of  winches. 

However,  a  more  reasonable  rate  of  1000  ft/min  pushes  the  strength  of  the  tether. 
Figure  U  shows  that  a  pay-out  rate  of  2000  ft/min  is  fast  enough  to  allow  the 
balloon  to  transverse  the  high  dynamic  pressure  region  before  it  is  blown 
down  range,  A  pay-out  rate  of  1000  ft/min  is  too  slow.  Figure  f>  further 
justifies  a  pay-out  rate  of  2000  ft/min,  A  1000  ft/min  pay-out  rate  is 
dangerously  close  to  the  ultimate  cable  strength,  while  a  2000  ft/min  pay-out 
rate  seems  to  maintain  a  factor  of  safety  of  1,6,  Figure  6  shows  that  the 
tether  is  in  no  danger  of  breaking  at  the  balloon  for  any  pay-out  rate. 

Figure  7  considers  the  dynamic  pressures  acting  on  the  balloon  as  it  ascends. 

Such  results  would  be  important  when  designing  the  balloon  structurally, 

A  graph  of  dynamic  pressure  versus  altitude  is  presented  for  each  run  in 

■the  results.  Their  importance  is  obvious  at  this  point  and  will  not  be  discussed 

further. 

The  second  set  of  Figures  (8,  9,  10,  ll)  compares  the  effects  of  three  different 
Summer  I  wind  profiles.  It  is  clear  (from  Figpre  8)  that  the  30  million  cubic 
foot  balloon  lacks  the  buoyant  lift  to  pull  it  through  the  high  dynamic  pressure 
region  in  a  95  percent  or  90  percent  wind.  However,  this  same  balloon  will 
ascend  in  a  75  percent  wind.  The  increasing  range  at  the  higi  altitudes 
for  the  75  percent  wind  can  be  attributed  to  the  increasing  tether  weight. 
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The  net  vertical  force  of  the  system  has  decreased  to  such  an  extent  that 
it  is  comparable  to  the  horizontal  drag  acting  on  the  system.  Consequently, 
the  horizontal  drag  force  effects  the  balloon's  motion  as  much  as  the  vertical 
forces.  Figure  9  also  demonstrates  that  a  balloon  of  this  size  will  break 
its  tether  if  launch  is  attempted  in  a  90  percent  wind.  But  an  ascent  in 
a  75  percent  wind  does  not  seriously  stress  the  tether.  Figure  10  indicates 
that  the  tether  strength  at  the  balloon  is  more  than  adequate  for  a  launch 
in  a  75  percent  wind. 

At  this  point,  it  was  reasoned  that  an  excess  of  helium  at  launch  would  give 
the  balloon  more  net  lift,  thereby  allowing  it  to  pass  through  the  high  dynamic 
region  before  it  was  blown  down  range.  The  results  of  a  Sumner  1-90  percent 
wind  profile  are  presented  in  Figures  12,  13,  lli,  l5.  As  was  expected, 
the  ascent  was  accomplished  as  can  be  seen  in  Figure  12.  However,  Figure  13 
clearly  shows  that  the  tether  undergoes  greater  stresses  and  actually  reaches 
its  ultimate  cable  strength  quicker  than  the  normal  balloon  does.  Figure  lli 
also  shows  a  large  increase  in  tension  at  the  balloon. 

The  next  endeavor  was  to  study  descent  trajectories  in  Summer  I  wind  profiles. 
First  consideration  was  given  to  variable  descent  rates  in  a  90  percent  wind 
(Figures  16,  17,  18,  19),  It  should  be  noted  that  the  balloon  used  in  this  run 
is  the  30  million  ft^  balloon  with  20  percent  of  the  helium  valved.  This  would 
result  in  a  buoyant  lift  which  wpuld  support  the  tether  payed-out  when  the 
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balloon  i.  at  5o,OCO  ft  in  Figure  8.  Figure  !6  si,ly  shows  that  .  >w 
w^oh  rats  allows  the  baUocn  to  be  blown  down  range  a  greater  distance 

than  a  fast  .inching  rate.  This  is  to  be  expected.  Figure  „  loaves  no 

ho*t  that  an  attested  descent  at  any  rate  would  sttess  tte  tether  beyond 
its  strength. 

figures  20,  21,  22,  23  compare  the  descent  of  the  normal  balloon  with  the 

one  discussed  above. (20  percent  helium  valved)  A  Cnmn  + 

varvedj.  A  computer  run  with  33.3  percent 

of  the  helium  valved  was  also  ran  and  plotted.  Tension 

P-LovwjQ.  iension  values  are  not  availabe 

for  this  run. 

In  order  to  demonstrate  a  successful  descwnt  trajectory,  .a  7S  percent  wind  was 
used.  Figures  21,,  25,  26,  2?  show  that  it  is  possible  to  retrieve  m  this 
if  a  slow  winching  rate  is  e^loyed.  Sven  at  200  ft/, in  the  factor  of 
safety  approaches  about  1.3  at  one  point  along  the  trajectory. 

A  second,  larger  balloon  wae  si*lated  for  winter  wind  conditions.  Observing 
that  a  Winter  1-25  percent  wind  profile  is  similar  to  a  S^r  I-7S  percent 

Ith  re;  11  13  r9aS0nable  *°  e3q39Ct  Sirailar  results-  Comparing  Figure  28 
gure  k,  both  trajectories  of  2000  ft/ndn  are  very  close.  Also,  Figures 

29  and  5  predict  about  the  same  factor  of  safetv  u 

s  fety  at  the  winch  throughout  the 

trajectory. 
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After  valving  out  19  percent  cf  the  helium,  a  descent  trajectory  was 
attempted  from  50,000  ft.  Figures  32,  33,  3U,  35  show  the  results.  As 
in  the  previous  ascent,  this  trajectory  compares  favorably  to  a  descent  in 
a  Sumner  1-75  percent  at  200  ft/min,  except  that  the  dynamic  pressure  is 
50  percent  higher. 
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The  results  of  this  preliminary  study  indicate  that  it  is  feasible  to  launch 
and  retrieve  high-altitude  balloons  under  certain  conditions.  It  has  been 
shown  that  a  30  million  ,f t3  balloon  can  be  launched  in  a  Sumner  1-75  percent 
wind  profile  at  2000  ft/min,  and  retrieved  in  the  same  wind  at  -200  ft/min 
with  20  percent  of  the  helium  valved  off.  Although  launch  and  retrieval 
are  possible  in  the  Winter  I  wind  profile,  it  is  more  restricted.  The 
79  minion  ft3  balloon  in  the  Winter  I  25  percent  wind  show  results  similar 
to  the  30  million  ft3  baHoon  in  the  Sumer  I  75  percent  wind. 

The  results  of  this  feasibility  study  can  be  used  as  preliminary  design 

criteria  for  components  of  the  tethered  baHoon-systems-banoon,  tether 
and  winch. 
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SECTION  IV  -  RECOMMENDATIONS 


Since  this  report  presents  the  results  of  a  feasibility  study,  the  only 
immediate  conclusion  reached  is  that  such  a  system  can  achieve  the  desired 
objective  of  deploying  tethered  balloons  through  high  wind  velocity  regions 
and  retrieving  them.  Since  the  physical  description  of  the  balloon  and 
tether  are  somewhat  idealized  (Appendix  B-2  and  B-6),  no  attempt  was  made  to 
study  the  system  in  detail.  It  is  recommended  that  a  balloon-tether  system 
be  designed  and  its  aerodynamic  properties  determined.  With  such  a  system, 
a  detailed  parametric  study  could  be  achieved. 

Further  recommendations  (relative  to  the  computer  simulation)  are  made  at 
the  end  of  Appendix  B. 
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Equations  of  Motion 


1.  General 

Using  Lagrangian  techniques  and  assuming  the  system  can  be  simulated 
by  a  tether  of  "N"  straight  links  with  a  balloon  hinged  to  the  top 
link,  two  differential  equations  can  be  derived.  The  position  of  the 
winch  on  a  flat  non-rotating  earth  is  considered  to  be  the  origin  of 
an  inertia  coordinate  system.  The  first  equation  expresses  the  motion 
of  the  balloon  about  its  hinge  point;  and  the  second  expresses  the 
motion  of  any  link  about  its  lower  connection  point.  The  derivation 
and  notation  follow  very  closely  the  work  of  Professors  Odin  R.  Elnan 
and  Carl  T.  Evert  (Reference  l). 

Lagrange's  differential  equations  of  motion  (Reference  2)  for  the 
system  are: 


where  the  generalized  coordinates  are  0  $nd  0  (r  -  1,  N),  the 

generalized  forces  are  FQ  and  (r  ■  1,  ...,  N)  and  T  is  the  total 
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kinetic  energy  of  the  system.  It  should  be  noted  that  the  coordinate  0 
could  be  chosen  as 

2.  Cable  Representation 

The  tether  (as  shown  in  Figure  A-l)  is  represented  by  "N*  links,  each 
of  constant  length.  However,  the  length  as  well  as  the  mass  of  each 
link  is  assumed  to  be  a  function  of  the  winching  rate  and  cable  profile. 
Each  link  is  considered  rigid,  and  each  hinge  fiictionless.  Therefore, 
all  forces  but  no  moments  can  be  transferred  through  a  hinge. 


Pn 


Figure  A-l  -  Tether  Model 
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The  angle  that  each  link  makes  with  the  horizontal  (measured  positive  for 
counterclockwise  rotations)  is  ^  (i  *  1,  N).  It  is  assumed  that  all 

external  forces  acting  on  a  link  can  be  concentrated  at  the  geometric 
center  of  the  link.  This  Introduces  an  obvious  error  for  both  the  aero¬ 
dynamic  and  gravitational  forces.  However,  since  the  ability  to  locate 
the  center  of  pressure  and  center  of  mass  for  each  link  contributes  to 
the  complexity  of  the  derivation  and  computer  simulation,  it  seems 
reasonable  to  ignore  the  error  at  this  point. 

A  typical  link  is  shown  in  Figure  A-2. 


Figure  A-2  -  A  Typical  Link 
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P^  denotes  the  mass  center  of  the  HiHth  link  and  is  assumed  to  be  a 
distance  of  i?c/2  from  the  hinge  point 


For  one  link: 


icx 


Pl*  ^C1  *  -*/  ^°1  — 

•  “2  el  T^el 


(3) 

(U) 


^Pl  -  ^cx  K  ®i  +  ^c-l  ®i 


(5) 


where  is  the  position  vector  from  PQ  to  P^ 


For  two  links: 


r2 


vP2  -vPl  +  v 


P2*/*l 


(6) 


p 

■2 


where  V  denotes  the  velocity  of  P^*  relative  to  Pi_1 

Vpi  i  -*/  ^c2  -s* 

-  F”  ^2  ®2  +  ”F  2 

A  / 

^  e^  +  /Cl  ^  +  “T^  K  ®2  +  “ F  ®2 


(7) 

(8) 


In  general: 


r  -  V 
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P  *  n-1  nZ^  .  Pc  .  /c 

-  £  4  4  «i'+  £  i0i  ri  *  V2  4  *n  *  -T  S 

i-1  i-1 


(9) 
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i-l 


i-1 


(10) 


Therefore,  the  kinetic  energy  of  the  tether  is: 


h 


kc  -  [t  “Cjj  (v  )2  +  ^  Ich 
n-1 


if] 


(11) 


If  each  link  is  assumed  to  be  a  uniform  thin  rod  the  moment  of  inertia  is: 

it 


E«n*%  U 


(12) 


T°  '  fml  ^  P”  )2<TC%^4 


(13) 


Other  useful  relations  to  be  used  later  are: 
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3.  Balloon  and  Harness  Representation 

The  balloon  harness  (denoted  by  ^  «  radius  of  balloon)  Is  considered  to 
be  normal  to  the  balloon's  longitudinal  axis  at  all  times.  The  balloon 
has  "directional  mass"  and  Mg  along  and  normal  to  the  harness  and 
concentrated  at  the  center  of  gravity  of  the  balloon  (Q).  Since  the  balloon 
is  assumed  to  be  an  expanding  sphere,  the  "directional  mass"  and  Mg 
are  equal  and  concentrated  at  the  geometric  center  of  the  sphere  as  shown 
in  Figure  A-3, 


«SF:  ENGINEEKING  MOCEDUtE  S-017 


GOODYEAR  AEROSPACE 

CORPOH4TIOM 

GER-2371U 


The  kinetic  energy  of  fee  harness  plua  balloon  is: 


0  2 


Vk-St-e***  ♦  4  ^  c*rti)2  .ihm  )2 ,  i  h  jf 


where  is  fee  velocity  of  Q 


balloon  about  Q. 


and  Ifi  is  the  mass  moment  of  inertia  of  the 


It  follows  that: 
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Using  equations  (?)  and  (22) 


p  #  n-1 

— **  n 

T 
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Using  equations  (17)  and  (22) 
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Using  equations  (l8)  and  (22) 
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Therefore  may  be  written  as: 

34 
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It  should  be  noted  that  the  time  derivative  of >^h  and  m^  has  been  excluded 
in  equation  (32).  Since  is  the  radius  of  the  balloon,  its  mass  is  zero 
and  therefore,  also  its  time  derivative.  The  radius  does  change,  however; 
but  its  time  derivative  is  small  when  compared  to  A  and  is  therefore  ignored. 


“he  time  derivative  of  equation  (31)  is: 
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Referring  to  Figure  A-3: 
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Expanding  all  the  terms  in  equation  (35)  results  In: 
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+  X“  0n  4  4  sin  <4  -  0r)  +  “r“  4  4  003  (4  “  0r) 

-  /cr  0r  cos(0n-0r)+  ^T“  ir  sin  (0n  “  0r)j 

r-l  jj  r-1  • 

+mcr  [  4  0i  4  sin  (0i  -  0r>+  %  K  — 5“  cos  (0±  -  0r) 

1-1  1-1 
r-l  r-l  • 

-«cr  {_  ici  -y1  0r  c°3  (01  “  0r)  +  %  ^ /ci  “T”  3in  "  4> 

1-1  1-1 

'  /N  N  \ 

+mh  -  ( £4  4  icr  0r  +  {  4  irj  sin  (4  -  4} 

-  >i-l  1-1  ' 


i  K  ir  -  ^  4  4  4  008  (4  -  4>  “  ®4  0r  003  (0r  +  e) 


-iA  *Jb  sin  (Q  +  0r) 
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+ml 


N 


“  cos  (4  +  e)icr  4  sin  (9+4) 


L  i-1 
N 


4) 


fM2 


-  £  ic±  sin  (0  +  4)  i£r  4  sin  (9  + 
i-1 
N 

+  {kk  COS  (4  +  0)  £cr  COS  (0  +  4) 
i-1 
N 

*{k  sin  (0  +  l±)Jcr  cos  (e  +4) 
i-1 
'  N 

/ 4  4  sln  (4  +  e)Ar  4  c°s  (0  +  4) 

i-1 
N 

‘  /  b±  cos  +  to£r  4  cos  (e  x  4^‘A  ®  ir  4  cos  (e  +  4) 
i-1 

N  .  N 

+  (L±  4  sin  C4  +  0)  4  si*  (0  +  4)-  /i  c°*  (e+4)i  sin (0+4) 

i-1  i-1  1 


/h  0  4  sin  (e  +  4) 


(36) 
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r  „ 

L  9 


r  8 

L  § 


Similarly, 


4^-*M  +VQ 

a  0  n  ^0  1  *  *N+1  ^  0  N+l 


_^N+1  ]' 

'  a  e  J 


Expanding  all  the  terms  In  equation  (37)  gives: 


-frm\ 


N  f)  "  ■  D 

-  /ic i  4-^f » °°»  <° +  4)  -  £  Zc±  hr 8  3in  (8  * 


♦M/4  4  cos  (j^  +  0)1  -  4  /Ci  0±  sin  (0  +  0i) 


;±  sin  (8  +  fUf  j  -A±  i>±  sin  (0  +  t>±) 


+  7  ^ci^iC08  (e  +  4^  ly  ioi  cos  (e  +  0±) 


+  /  /  sin  (0  +  0a 


.1/9  3 


Ci  COS  (P  •  0±) 


*  f  k  0j_  sin  (0j_  +  9)  /  4/ci  0i  cos  (0  +  0j_) 


'  /  /  li-1 


(continued  on  next  page) 
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Substituting  equations  (33)  and  (38)  into  equation  (l),  cancelling  like 

t  • 

terms  and  solving  for  6  yields: 

f  /  N  N 

6  ‘  j  Fe  "  0  *b  ♦  |j%  *T~  l  003  (4  ♦  e)  +  £  4  sin  (4  *  ejj 


*  2  ^0  "hih  ^  4  sin  (e  *  0.)  +  «2  -  M,]. 

1*1 

‘  "  /  »  N 

/0  2  4  Sin  (8  ♦  4)  (ij.  2  4  00s  (4  <•  8)  +j[  ^  ain  (0  +  0±)  j 
W  1-1  i-1  / 


(continued  on  next  page) 
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Substituting  equations  (32)  and  (36)  into  equation  (2),  cancelling  like 
terms  and  solving  for  $r  yields: 


K 


cos  (0±  -  0r) 


n-1  n-1 

+  k^fclc  £  sin  (0±  -  0r)  -  ^  sin  (0±  -  0J 

i-1  jl-3. 


N  n-1 

'1*  [  %  {  'k  008  wi  -  V 

n-r+2  i-r+1 


i-1 


cos  (0±  -  0r ) 


(continued  on  next  page) 


tEF:  ENGINEERING  MIOCECUIE  S  0I7 


GOODYEAR  AEROSPACE 

CORPORATION 

_ GER-1371U 


involving  are  set  equal  to  zero. 


5.  The  Generalized  Forces 

The  generalized  force  F0  is  given  by  ; 


P*  N 

-  N+1  _ 

_ ,  s'  3v  h 

9  "»♦!  Je  «'tr" 


F  -  F 


Jv 


p  * 

4  i-Jtl  ♦  ■ 

C  n  ^  o 


n«l 


where 


Oa) 


F0  is  the  external  torque  acting  on  the  balloon 

_ k 

FN+1  is  aPPlied  force  acting  on  the  harness 

Fq  is  the  applied  force  acting  on  the  balloon 

Fn  is  the  applied  force  acting  on  the  link 

M  is  the  applied  moment  about  the  center  of  mass  of  the  balloon 

-Pn* 

From  equation  (9),  V  is  not  a  function  of  0. 


The  applied  moment  about  the  center  of  mass  is: 


H 


)  e 


6  +  M 
s 


(12) 


where ; 

M  is  the  aerodynamic  moment 

3 

the  aerodynamic  damping  moment 
static  moment  due  to  buoyant  lift 

A  -23 


Mg  is  the 


Since  the  balloon  is  spherical,  the  buoyant  force  is  concentrated  at 
the  geometric  center,  and  Mg  ■  0.  A  spherical  body  also  has  a  zero 
aerodynamic  moment*  If  viscous  effects  are  ignored,  the  aerodynamic 
damping  moment  is  zero*  Therefore,  equation  (1*1)  simplifies  to: 


-»F* 

-  N+i 

F  - F  . 

0  N+l 


+ 

it  q  Js 


Fq  -  HH1+BH2 


(U3) 


Fq  represents  the  force  vector  acting  on  the  center  of  gravity  of  the 
balloon  and  is  given  by: 


m 


where  H  is  the  horizontal  force  and  B  is  the  vertical  force  as  shown  in 
Figure  A-l*. 

Vg  ■  Balloon  velocity 
V  -  Wind  velocity 
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The  forces  that  contribute  tc  B  and  H  are  the  buoyant  force,  the  weight 

of  the  balloon  struck,  the  weight  of  the  helium  Inside  the  balloon  and 
the  drag  forces. 

The  total  drag  force  is  Dg  -  qg  Sg 

cDg  is  assumed  to  be  constant  and  equal  to  .15 
SB  varies  and  is  equal  to  -^2  ^ft2 

%  is  the  dynamo  pressure  acting  on  the  balloon  relative  to  the 
air  xu.  lb/ft2 


In  this  case,  not  on*  is  the  balloon  nortl*  but  a  variable  wind  is  acti* 
on  it.  The  wind  velocity  is  a  function  of  altitude.  At  ary  tine,  the 
wind  velocity  is  considered  constant  and  equal  to  the  velocity  of  the 
wind  at  the  altitude  of  the  center  of  gravity  of  the  balloon.  This  assertion 
eliminates  an  aerodynamic  moment  which  probably  does  slightly  effect  tte 
notion  of  large  balloons.  If  the  balloon  is  moving  vertical*  witt  a  velocity 

a  and  horioontal*  with  a  velocity  i  am  the  horiaontal  wind  velocity  is 
Vg  than  the  dynamic  pressure  is: 


SB 


■  i  f  -  V2  ♦  a2] 


0*6) 

The  drag  can  now  be  brohen  into  two  components,  one  in  the  direction  am 
one  into  the  N2  direction. 


A-25 


V — 


I 

r 

r 

r 

r 

i 

c 

r 

r 

E 

e 

E 


I  ?r 

r  s 

L-  i 

,  *  -O  is 


3  z 

2  » 
k  2 
-1 
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Referring  to  Figure  A-l*: 

*  -  tan"1  [z/  (x  -  VgVj 
Then  the  forces  on  the  balloon  ares 


B  -  lg  -  %  -  Wgg  -  %] 


B 


H  - 


where: 


lB 


1^  is  the  buoyant  force  of  the  displaced  air 

Wjj  is  the  weight  of  the  helium 

Wg  is  the  weight  of  the  balloon  structure 

is  the  drag  force  in  the  vertical  direction 
(positive  down)  and  is  equal  to: 


\  “  %  %  sin* 

Du  is  the  drag  force  in  the  horizontal  dire  tion 

nB 

(positive  to  the  left)  and  is  equal  to: 

X  a  %  %  cDb  C0Se( 


Since  the  harness  has  no  weight  or  aerodynamic  forces  acting  on  it, 


F  -  0 
N+l 


(U7) 


(lifl) 

(k9) 


(50) 


(51) 


(52) 
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Therefore,  when  substituting  equation.  (1J),  (W)  ^  (52)  tato  (lt3) 
the  generalized  force  on  the  balloon  becomes: 

Fe  -  (H  Nx  ♦  B  »2)  (J  ^  , 


(53) 


Simplifying 


B  sin  0  +  H  cos  ej 


m 


j 


From  Figure  A-l*,  it  is  obvious  that  FQ  is  the  external  torque  about  point 
PN.  (positive  torque  clockwise). 


J 


O 

U  * 

*  o 
•  z 


^  I 

>  <  o 

*: 

2  X 


The  generalized  force  F^  is  given  by: 

*r 
N 

,^PN+1  _Q 

F* 


K 


P*  „* 

/  —  ,-SPi  ^  — *  PN+1 

7  f  .11 —  +  f  .ii _  , 

^  1  JKr  N+1  ^T~ 

i“l  r 


Q'  aJL 


(55) 


N 


only  unevaluated  terms  are  £  f  uhich  are  the  applied  force. 


acting  on/c^  shown  in  Figure  A-2. 


i-1 


_  P *  p1 

Fi  =  F  1  +g‘  1 

T5  # 


(56) 


where  G  is  the  gravitational  force  of 


Pi* 


-  «ti  8  »2 


(57) 


li 

u 


A-27 


£ 


and  F  is  the  normal  drag  force  on 
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-  \  sln2  4  -  DJCicos2  4  4' 


.2  4  ■*' 


As  with  the  balloon,  each  link  is  moving  and  a  variable  wind  is  acting  on 

it.  At  any  time,  the  wind  velocity  across  a  link  is  considered  constant  and 

equal  to  the  velocity  of  the  wind  at  the  altitude  of  the  center  of  mass  of 

the  link.  It  should  be  remembered  that  •‘he  center  of  mass  was  assumed  to 

be  at  the  geometric  center  of  each  link.  If  the  center  of  mass  of  a  UnV 

is  moving  vertically  with  a  velocity  z  and  horizontally  with  a  velocity  x_ 

ci  ci 

and  the  wind  velocity  at  the  center  of  mass  is  V„  ,  then  the  drag  components 


\  -  *  Vi  sCi  K  - 

"  2  Cdc  'h  S=1  [\] 


The  factor  sin^  0^  and  cosc  0i  in  equation  ($8)  take  into  account  the 

projected  area  in  the  horizontal  and  vertical  directions  and  then  rotate  the 

horizontal  and  vertical  drag  components  into  the  direction  normal  to  A  . 

i 

Cq  is  assumed  to  be  constant  and  equal  to  1,2 
c 

Sc^  is  equal  to  the  average  diameter  of  J[c ^  times  Jc 
is  measured  at  the  center  of  mass  of  Jc 
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Therefore,  when  substituting  equations  .(lU,  (I4H),  (52)  and  (56)  into 
equation  (55),  the  generalized  force  on  each  link  becomes: 


(61) 


+  (H  Nx  +  B  N 2)Jc  £ 

Simplifying: 

\  '  '  4-  V  ^  4  -  4-  X  oos2  -  %  E  4-  cos  4 

r  p  A  t  r 

N 

*L  -  %0±  ^n2  4  oos  (4  -  4)  .  cos2  4  cos  (4  -  4) 


i-r+1 


-  raci  g  cos  0rj  +/c  [*-  H  sin  ^  +  B  cos  ^ 


(62) 


From  Figures  A-l,  A-2,  and  A-l*,  it  can  be  seen  that  is  the  external 
torque  acting  about  (positive  counterclockwise).  * 
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Computer  Program 


1.  Wind  Profile 

The  wind  profile  for  any  simulation  is  considered  to  be  a  function  of 
altitude  as  shown  in  Figures  1  and  2.  Any  one  of  these  curves  can  be 
approximated  by  twenty-four  points  and  twenty- three  line  segments.  The 
twenty-four  points  are  stored  in  an  array  in  the  computer  and  the  program 
interpolates  between  any  two  points.  Since  the  wind  profile  curve  is  well 
behaved,  5,000  ft  increments  seem  to  be  more  than  adequate.  However,  any 
altitude  increment  may  be  used,  and  the  increments  can  be  of  varying  size. 

2.  Tether  Profile 

Tether  profiles  used  in  this  simulation  are  shown  in  Figure  B-l.  Both  the 
diameter  and  weight  per  unit  length  can  be  approximated  as  a  series  of  linear 
functions  of  length,  and  is  therefore  represented  by  thirty-two  points  and 
thirty-one  line  segments.  The  computer  program  interpolates  to  find  the 
diameter  and  weight  per  unit  length  at  each  end  of  a  link  and  then  considers 
each  link  to  be  a  right  circular  cylinder  with  diameter  and  weight  found 
by  averaging  the  end  conditions.  Obviously,  the  greater  number  of  links  used 
the  more  accurate  the  approximation.  It  has  been  calculated  that  a  three 
link  representation  of  135*000  ft  of  tether  overestimates  the  tether 
weight  by  less  than  3.5^. 
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3.  Winching  Rate 

At  any  time,  the  winching  rate  is  considered  constant.  However,  the 
program  will  consider  discontinuous  changes  in  the  rate  as  a  function  of 
time  as  shown  in  Figure  B-2. 
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it.  Translational  Dynamics  of  Balloon  and  Tether 

Although  angular  coordinates  represent  ^  be£t  of  3oMng  ^ 

equations  of  motion,  Unear  coordinates  must  he  calculated  to  find  the 
forces  acting  on  the  balloon  and  links. 


a.  Balloon  Itynamics 


The  following  equations  represent  the  linear  displacewnts, 
velocities,  and  accelerations  of  the  balloon: 


L{  cos  4  4  si 


sin  6 


L  ^  sin  0±  +Jh 


cos  Q 


1c  ^  cos  fa.  -  jjc  ^  sin  0±  ff± 


h  cos  e  e 


«  N 

2  -  ic  ^  sin  4  ♦  /,  ^  COS  4  4 


sin  9  9 


I 

i 

I 


i 

I 

t 


r 

r 

r 
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z  «  2 


-  y  N  N 

4  £  cos  4  4  -4  /  sin  4  42  +/c  /  ooa  4  4 
1-1  i=l  ^ 


-/  cos  0  02  -/ 


sin  0  0 


(67) 


x  *  -2 


.  *  N  N 

4  £  sin  4  4  -4  /  cos  4  42  .  /c  /  sln  4  4 
W  i-1  f-i 


-/h  sin  0  02  +J?h  cos  0  0 


(68) 


b.  Tether  Dynamics 

The  following  equations  represent  the  LLnSar  vertical  displacement, 
velocities  and  accelerations  of  the  geometric  center  of  the  "r"th  link. 


r-1 


z°r  "  A  /  sin  4  +  —■  sin  K 


(6?) 


i-1 

r-1 


r-1 


cos  4  -  A  {  sin  0^  0^ 


i-1 


i-1 


L  J  (j 

~2~  cos  $r  "  -I2  sin  &  4 


(70) 
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r-x  r-1 

acr  *  Iz  ^  sin  4  +  A  ^  cos  4 
i“l  i-l 

i  0 

+  sin  0T  +  -|£  cos  4  ^ 


cr  ' 


/  *  *  r-i 

£  cos  *  *  4Z=  /  -  4  4  -Zc  /  ■*.  a  V 


X  cos  4  4  +  cos  4  4  -  •— 


r  T  sin  4  42 


A  r-1  r  ^ 

%  -  sin  4  4  -4  £  cos  4  42  ./c  /  sin  ^  jj 

i-l  *  ml  L- 

11  i-l 


-ic  sin  44-^ cos  4  42  -  A 


sin  4  4 


5.  Tension 


In  an  effort  to  show  ho»  «ch  loading  the  tether  must  be  able  to  accommodate, 
the  tension  at  each  hinge  point  was  calculated.  Referring  to  Figure  B-3,  the 
tension  at  the  balloon  hinge  point  is: 

£P*  *  V*  ■  B  -  a  -  Fy 


4jj  *  B  -  PL  -  4  z' 


B-6 
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Ux  "  ^  "  H  "  % 

%  "  H-l^x 


(V  +  <V* 


where : 


31X1  /FX  are  the  sum  of  the  vertical  and  horizontal 
forces  respectively  acting  on  the  balloon. 

and  are  the  vertical  and  horizontal  components 
of  tension  on  the  N**1  link. 


PL  is  the  payload  weight  in  lbs. 


Tn  is  the  tension  in  the  tether  at  the  balloon  hinge. 


B-7 
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Figure  B-3  -  Horizontal  and  Vertical  Forces  Acting  on  Balloon 


Figure  B— 14.  shows  the  "N"4«h  llnk  and  all  the  forces  acting  Qn 


■  %  S  *  \  -X'X- 

N-l 

(79) 

Fv 

N-l 

■  x  -  \  -  x  -  «oN  i;N 

(80) 

•  * 

■  x  »  Fu„  +  n  «.  f 

CN  HN  % 

(81) 

%-i 

“  Fw  +  °h  -  x’ 

%  hn  cn  °n 

(82) 
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Vl  “d  %_!  3re  4116  ,ertical  ^  horizontal  components  or 
tension  at  the  top  of  the  "N-l"  link. 

WTN  is  the  weight  of  the  *'N*,th  link  in  lbs. 

In  general  for  the  ''r"th  link 


The  tension  at  the  winch  is  found  by  setting  r-1 

\  \  +  \  ^  xCl 

Fvo  "  \  '’\-Wri-mc1 

%  ■  i  <v2  +  v 


m 

(85) 

(86) 


(67) 

(88) 

(89) 


6.  Physical  Aspects  of  the  Balloon 

As  mentioned  previously,  the  balloon  is  considered  to  be  a  thin  spherical 
shell  expanding  as  it  ascends.  Two  basic  balloon  sizes  were  used.  The 
rirst  one  emended  to  30  x  vfi  «5  at  100,W0  ft  altitude,  and  the  second, 
to  79  x  10  ft3  at  the  same  altitude.  The  structural  weichts  of  the  balloons 


B-10 


GOODYEAR  AEROSPACE 

CORPORATION 

GER-I371U 


are  8600  lbs  and  2b,600  lbs  respectively.  The  volume  at  any  altitude 
is  equal  to  the  volume  at  sea  level  divided  by  the  atmospheric  density  ratio. 
The  weight  of  helium  (96.5/6  pure)  in  the  balloon  is  equal  to  the  weight 
density  of  helium  times  the  volume  of  the  balloon.  For  the  30  million  ft3 
balloon,  the  weight  is : 

\  -  (.01297)  (U19000)  =  5h30  lbs  (90) 

-3 

For  the  79  million  ft  balloon,  the  weight  is: 

WR  -  (.01297)  (1,102,000)  «  lh,300  lbs.  (92.) 

The  payload  is  assumed  to  be  attached  to  the  balloon  at  the  tether 

hinge  point.  The  moment  of  inertia  about  the  c.g.  of  the  balloon  takes 

into  account  the  payload  but  not  the  helium  inside  the  balloon  or  any  of 

the  air  outside  the  balloon.  This  assumes  that  the  shear  forces  are  negligible . 

IB  •  (2/3  'Ju3  *  P L)Jhye  (92) 

The  directional  masses  of  the  balloon  (Mj_,  M2)  are 

h  m  h  “  (WH  +  \  +  PL)/g  +  IPo  Vo  (93) 

where : 

\  (  Vq  in  the  apparent  mass  of  the  balloons 
f  is  the  atmospheric  density  a  t  sea  level 

V  is  the  balloon  volume  at  sea  level 
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and  are  equal  because  the  balloon  is  spherical. 


7.  Numerical  Integration  Technique 

The  two  dynamical  equations  that  must  be  solved  are  (39)  and  (Uo).  These 
equations  can  be  written  in  the  form: 


•  •  ,  |  |  ••  »»  i*  *  j  *  • 

0r  *  6  9*  $r-l*  ^r+l* ###*  $!*•••*  *k) 


(9U) 


It  is  clearly  seen  that  equations  (91*)  are  strongly  coupled.  To  reduce  the 
complexity  of  the  numerical  integration,  all  coupling  terms  of  second  order 
are  eliminated.  Preliminary  computer  runs  show  that  angular  accelerations  are 
one  order  of  magnitude  less  than  angular  velocities.  Therefore,  equations  (91;) 
reduce  to  the  following; 


it  •  * 

0  *  f]_  (®j  $]_»•••*  ^  ) 

Equations  (95)  are  solved  by  Runge-Kutta  (l;-^  order  accuracy)  numerical 
integration  techniques.  Due  to  the  long  running  times,  it  is  advantageous 
to  use  a  large  time  increment.  It  has  been  found  that  an  increase  in 
the  degrees  of  freedom  requires  a  reduction  in  the  time  increment.  To  help 
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alleviate  this  problem*  the  balloon* s  pitch  angle  is  held  constant. 

This  affect  on  the  dynamics  of  the  tether  was  hardly  noticeable.  No 
graphs  were  constructed  to  show  this  because  the  differences  were  in  the 
fourth  significant  figure. 
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8.  Deternination  of  Number  of  links 

As  stated  previously,  an  increase  in  the  degrees  of  freedom  decreases  the 
magnitude  of  the  time  increment  which  will  allow  a  stable  simulation.  Also 
the  computer  time  increases  linearly  with  an  increase  in  the  number  of  links 
simply  because  each  extra  link  involves  one  more  equation  of  the  form  of 
equation  (UO)  which  must  be  integrated.  Therefore,  an  effort  was  made  to 
determine  the  least  number  of  links  which  would  present  a  gocd  simulation. 

To  this  end,  a  tether  was  payed-out  in  a  summer  1-9 ($  wind.  A  winching 
rate  of  3000  ft/min  was  used  for  1500  sec.  This  allowed  the  balloon  to  pass 
through  the  high  dynamic  pressure  region.  Figure  B-5  shows  the  altitude 
versus  range  of  the  balloon  for  both  three  links  and  five  links.  Figure  B-6 
shows  the  tether  profile  at  500  sec,  1000  sec,  and  l500  sec.  These  comparisons 
justify  the  use  of  three  links  in  all  future  runs  unless  a  wind  profile 
containing  severe  cross  winds  is  used.  The  conclusions  that  three  links  is 
adequate  for  the  simulation  was  also  reached  by  Professors  Elnan  and  Evert 
in  Reference  1. 
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The  format  for  all  inputs  is  F10.0.  The  following  is  a  list  of  all 
needed  variables  with  appropriate  units: 


a.  Stop  Altitude  -  The  program  may  be  stopped  when  the  balloon  reaches 

any  given  altitude  (ft)  and  a  new  simulation  started.  This 
variable  is  good  only  for  ascent  simulation. 

b.  Stop  Time  -  The  program  may  be  stopped  at  any  given  time  (sec)  during 

the  simulation  and  new  simulation  began. 

c.  Stop  Tether  Length  -  The  program  may  be  stopped  when  a  certain  length 

of  cable  (ft)  has  been  winched  out. 

o 

d.  Acceleration  of  Gravity  at  Sea  Level  -  Ft/sec 

e.  Radius  of  Earth  -  Ft 

f.  Integration  Step  Size  -  The  program  is  started  with  a  relatively  small 

time  increment  (.5  sec)  to  allow  for  large  acceleration  to  damp 


g.  Printout  Time  -  In  the  majority  of  simulations,  it  is  not  necessary  to 

printout  after  every  integration.  Since  printing  time  wastes 
computer  time,  it  is  most  desirable  to  make  several  integrations 
before  printing  out  data.  This  input  variable  should  be  equal 
to  the  number  of  integrations  made  before  each  printout. 

h.  Number  of  Links 
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i.  Weight  of  Payload  Attached  at  Balloon  -  Tether  Junction  *-  Bb 

j.  Time  -  For  an  ascent  simulation,  the  time  is  arbitrarily  set  at  k  sec. 

The  program  then  starts  with  a  relatively  small  length  of  tether 
winched  out.  If  the  length  of  tether  were  zero  initially, 
oscillations  of  a  very  high  frequency  would  take  place  for  a  short 
period  of  time  because  the  moment  arm  is  very  small.  For  descent 
simulations,  the  time  (sec)  must  be  equal  to  a  value  which,  when 
multiplied  by  a  positive  winching  rate,  will  equal  the  total 
length  of  tether  payed-out  at  the  start  of  the  descent. 

k.  Initial  Tether  Length  -  Ft 

l.  Second  Integration  Step  Size  -  After  the  first  several  seconds  (sec  f) 

the  time  increment  is  increased  in  order  to  speed  up  the  calculations. 

m.  Second  Time  -  This  variable  is  set  equal  to  the  time  (see  j )  plus 

several  seconds  (10-20).  After  the  3econd  time  is  reached,  a  second 
larger  integration  step  size  is  used  (see  l)  because  accelerations 
are  relatively  low. . 

n.  Initial  Altitude  -  Ft 

o.  Drag  Coefficient  of  Cable  -  A  value  of  1.2  (cylinder)  was  used  for  this 

simulation. 

p.  Drap  coefficient  of  balloon  -  A  value  of  .15  (sphere)  was  used  for  this 

simulation. 

q.  Reference  Volume  of  Balloon  -  The  program  uses  the  volume  (ft^)  of  the 

balloon  at  sea  level  as  a  reference. 
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r.  Harness  Mass  -  Slugs  (zero  in  these  runs) 

s.  Balloon  Structural  Weight  -  Lb 

t.  Weight  Density  of  Helitun  at  Sea  Level  -  Lb/ft3 

u.  Initial  Pitch  Angle  of  Balloon  -  Rad  (zero  in  these  runs) 

v.  Initial  Pitch  Rate  of  Balloon  -  Rad/sec  (zero  in  these  runs) 

aa.  An  array  of  N  variables  representing  the  initial  tine  rates  of  change 
of  the  angles  of  the  N  links  -  positive  counterclockwise  -  Rad/sec 

bb.  An  array  of  N  variables  representing  the  angles  of  N  links  -  positive 
counterclockwise  from  the  right  horizontal  -  Rad 

cc.  An  array  of  eight  variables  representing  times  at  which  the  winching 
rate  is  changed  -  sec 

dd.  An  array  of  eight  variables  representing  winching  rates  corresponding 
to  (cc).  The  first  value  is  the  initial  winching  rate  and  corresponds 
to  zero  time.  The  second  value  is  also  the  initial  winching  rate  and 
corresponds  to  the  time  at  which  the  winching  rate  is  to  be  changed. 

The  third  value  is  the  second  winching  rate  and  corresponds  to  the  time 
at  which  the  winching  rate  is  to  be  changed  again,  and  so  forth.  -Ft/sec. 

ee.  An  array  of  twenty-four  variables  representing  altitudes  -  Ft. 

ff.  An  array  of  twenty-four  variables  representing  wind  velocities  corresponding 
to  the  altitudes  given  in  (ee)  -  Ft/sec. 

gg.  An  array  of  thirty-t»„  variables  representing  lengths  of  the  tether 
starting  at  the  balloon  -  ft. 
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hh.  An  array  of  thirty-two  variables  representing  diameters  of  the  tether 
corresponding  to  (gg).-  Ft. 

ii.  An  array  of  thirty- two  variables  representing  weights  per  unit 
length  of  the  tether  corresponding  to  (gg)  -  Ib/ft. 

10.  Outputs 

At  predetermined  time  intervals  (see  input  g.  )  the  following  data  is 
outputted  in  format  F9.2. 

a.  The  angle  formed  by  each  link  from  the  horizon  -  Deg 

b.  The  angular  velocity  of  each  link  -  Deg/sec 

c.  The  tension  at  the  top  of  each  link  -  Ih 

d.  Time  -  Sec 

e.  Altitude  of  Balloon  -  Ft 

f.  Range  of  Balloon  -  Ft 

g.  Horizontal  and  vertical  velocities  of  balloon  -  Ft/sec 

h.  Pitch  angle  of  balloon  -  Deg 

i.  Pitch  rate  of  balloon  -  Deg/sec 
Length  of  one  link  -  Ft 

k.  Winching  rate  -  ft/min 

l.  Wind  velocity  at  balloon's  altitude  -  Ft/sec 

m.  Harness  length  (balloon  radius)  -  Ft 

n.  Tension  at  winch  -  Lb 
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0.  Total  weight  of  cable  winched  out  -  lb 

p.  Horizontal  drag  force  on  balloon  -  lb 

q.  Vertical  drag  force  on  balloon  -  lb 

r.  Total  vertical  force  acting  on  balloon,  including  the  buoyant  lift, 
structural  weight,  weight  of  helium  and  vertical  drag  force.  This 
force  must  be  greater  than  the  weight  of  the  tether  plus  the  payload 
in  order  to  have  an  ascending  balloon  -  Lb. 

s.  Dynamic  pressure  acting  on  the  balloon  relative  to  the  free  stream 

2 

velocity  -  Lb/ft 
11.  Fortran  IV  Program 

A  listing  of  the  program  is  contained  in  APPENDIX  C.  The  functions  of 
the  main  program  and  the  subroutines  are  as  follows: 

a.  Main  Program 

1)  Define  the  two  functions  to  be  solved  (equations  (39)  and.(UO)). 

2)  Read  input  variables 

3)  Initialize  program 

U)  Check  if  end  constraints  (balloon  altitude,  time,  and 

tether  length)  are  met.  If  one  condition  is  met,  print  output, 
and  start  new  program;  otherwise,  continue. 

5)  Check  if  printout  time  is  reached.  If  it  is  reached,  calculate 
the  tensions  at  the  top  of  each  link  and  at  the  winch,  print 
output  and  continue;  if  not  reached,  continue.. 
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6)  Advance  the  generalized  coordinates  and  velocities  through  one 
time  increment  by  Runge-Kutta  integration, 

7)  Calculate  altitude,  range,  horizontal  and  vertical  velocities 
and  accelerations  of  balloon, 

8)  Calculate  time  rate  of  change  of  mass  of  each  link  and  time 
rate  of  ch singe  of  moment  of  inertia  of  balloon, 

b.  Subroutine  Head  (N) 

*  * 

This  subroutine  writes  out  the  headings  0^,  0^,  0^  0^ 

and  ThNS^ ,  •  •  * ,  TENS-^2  * 

c ,  Subroutine  Subrb 

1)  Find  the  atmospheric  density  surrounding  the  balloon  by  using 
an  "in-house"  subroutine  based  on  1962  US  Standard  Atmosphere. 

2)  Calculate  volume,  moment  of  inertia  about  c.g,,  radius  and  aerodynamic 
reference  area  of  balloon. 

3)  Calculate  horizontal  wind  velocity  at  balloon's  altitude, 

U)  Calculate  horizontal  and  vertical  drags  on  balloon  due  to 

relative  velocity  between  balloon  and  wind, 

5)  Calculate  directional  masses  of  balloon, 

6)  Calculate  static  lift  of  balloon, 

7)  Calculate  generalized  force  acting  on  balloon, 
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d.  Subroutine  Subrc 

1)  Identify  several  variables  with  the  appropriate  link. 

2)  Calculate  the  generalized  force  acting  on  the  “r"  link. 

e.  Subroutine  Suirfc 

l)  Evaluate  all  of  the  summations  found  in  the  ©equation  (39). 

f.  Subroutine  Surac 

1)  Evaluate  the  remaining  summations  to  be  used  in  the  &  equation  (bO). 

g.  Subroutine  Tether 

1)  Calculate  the  length  and  time  rate  of  change  of  the  links. 

2)  Calculate  the  mass  and  aerodynamic  reference  area  of  each  link. 

3)  Calculate  altitude,  horizontal  and  vertical  velocities  and 
accelerations  of  the  geometric  center  of  each  link. 

h)  Find  the  atmospheric  density  at  the  geometric  center  of  each  link. 

5)  Calculate  the  wind  velocity  at  the  geometric  center  of  each  link. 

6)  Calculate  the  horizontal  and  vertical  drag  of  each  link, 

12 .  Recommendations 

Throughout  APPENDICES  A  and  B  actions  and  approximation.,  were  made  so 
that  the  complexity  of  the  equations  of  motion  and  the  computer  program  would 
he  held  to  a  minimum.  Since  the  program  called  for  a  feasibility  study  of 
a  balloon- tether  system,  the  results  obtained  are  considered  meani-gful  and 
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fairly  accurate.  However,  if  a  complete  parametric  study  is  attempted, 
several  refinements  could  be  made. 

a.  The  location  of  the  center  of  pressure  and  center  of  gravity  of 

each  link  should  be  calculated  and  not  assumed  to  be  at  the  gt  metric 
center. 

b.  A  better  description  of  the  physical  and  aerodynamic  properties 
of  the  balloon  and  harness  should  be  defined. 

c.  A  variable  wind  profile  over  each  link  and  a  link  that  varies  in 
diameter  should  be  considered  in  the  computer  program. 

d.  The  weight  of  each  link  should  be  found  by  integrating  a  polynomial 
which  represents  the  weight  per  unit  length  rather  than  averaging 
end  conditions. 
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Program  listing 

DISK  OPERATn^  SYSTEM/366  FORTRAN  36CN-FQ-451  22 

DOUBLE  PRECISION  T,DT  . 

_  DIMENSION  AA(4»12)»BB(4) » PH  IDEGI 12)  *P  HIDDEI  12)»MC1(12)»DDD(12) 
DIMENSION  2C(12) , SC( 12 ) t RH0C1 12) . VCI 12)»VGC( 12), XCD( 12) ,ZCD( 12) 

__  DIMENSION  VF!12)fHF.!12) _ 

REAL  MC{12)tMCD(12)tLHtLC.MlfM2,MHtLCD.IB.LCLC»MCRD,MCRtLCLH, _ 

1LC LCD  t LLCD ( 8)  « LL  ( 32) ,MM1 1 8 ) , MM2 (8 ) t LCSL tNCR.NC ( 12) , I BD, I B1,MC1, LS 

2,NCV(12)»NCVR  _  .  _ 

COMMONT»  Sl»S2»S3»S4fS5»S6» S7, S8 , S9t S10.S11 , S12 , S13 ,S14 » S15 »S16 ,S17 
ltS18tS19jS20»S21tS22»S23»S24»SSl.SS2»SS3»SS4»SS5»PHID(12) »PHI I 12 ) » 

_ 2  THE  , PH  I OD ( 1 2 ) , PH  I R , $ B ,  DYPRB  »  C.QB*  Hi.BxVQL  B  ,.R  HOj  G  «  W  3 » CM . CMQ*  CAS*  V(L* _ 

3VB , THEDfFTHE  »FPHI ,AN » AR , I R , X , Z , XD , ZD , N»TT ( 8 ) t RREN ( 8 ) , PHI RD, VOLO. 
4CCD(8)«VVG(24).ZZ( 24  ) , WDHSL. WBS. DD( 32) ,NC» SC.VGC.VC .NCV.NC VR . 
5MC.MCD.LH.LC.M1 , M2 , MH , LCD , I B , LCLC , MCRD, MCR, LCLH. LCLCD. WWPUL ( 32  ) . . 
6LLCDtLL»MMl.MM2.TETH,S25fS26,WHtWT(12 ) , NCR, D ( 1 2 ) » CDC» WTCT , TENS( 12) 
_7»TENSW.DVB»PL«ZCOD( 12) ,DVC , XCDDI 1 2 ) 

C _  _  _ _ _  _  _ 

C  EQUATIONS  OF  MOTION 

c 

Jioob  FTHEDD(THED)=(FTHE+(M2+MH/2.)*(S1+S2)*LCLH 

l  +  ( 2.  *M2+MH ) *LH*LC  0*S3+( M2-M1 ) *{ LC*S3* ( LC*S4+LCD*S5 )-LCD*S6* (LC*S4+ 
2LCD*S5) )-THED*IBD)/( I3+LH*LH*( MH/3.+M2) ) 

2000  FPHIDD(PHID)  =  (FPHI+LCLC*(-SS1+SS4.  +lrS17-S19+S?l-MCRD*S9-Mj;R*Sl_0 
l+MCR*S12)/2.-S25*S10-SS5-MH*Si6+MH*S22+Ml*THED*SIN( THE+PHIR)* 

2S4— M 1*C0S  ( THE+PH I R )*( S23+S1 5-S7-THED*S3 )-M2*CGS ( THE+PHI R) *THED* 

3S3— M2*  SINITHE+PHIR )*( S 13+S 14+S 1+THED*S4)-MH*  S10-MCRD*PHIRD/3 • )  + 
4LCLC0*(-2.*SS2-SS3-S18-S20/2.-MCR*S9-MCRD*Sll/2.-2.*MCR*PHlRD/3. 

5-2. *MH*S8-2.*M1*C0S ( THE+PHI R)*S4+M1*S INI THE+PHI R )*THED*S5- 

_ 6M1*C0S(  THE  +  PHI  R)*THED*36-M2*SJNJIHE  +  PJHIRJ*(2.*S3+THED*S5J± _ 

7M2*CCS (THE+PHIR )*THEC*S6 ) 

8+LCLH*THED*THED*C0S (THE+PHI R)*(MH/2.+M2) ) / (L CLC* ( S25+MCR / 3. + 
9MH+M1*C0S (THE+PHI R)*C0S (THE+PHIR ) +M2*SIN( THE+PH I R ) *SI N ( THE+PHIR) ) ) 

50  FORMAT ( 14X ,4HT IME.6X, 3HALT. 7X.  5HR ANGE , 5X, 6HH0R VEL , 4X, 7H VER TVEL . 3X » 

_ 15HTHETA.5X.3HTHE-RATF,  2X, 4HL INK , 6X, 7HP0-R ATE , 3X , 7H WIND VE L , 3X, 2HLH , 

_ 2  8X,5HTENSW/14X,7HWTCABLE,3Xt5H0RAGHj5Xj5HDRA')V,5X,8HVFRTF0RCJLZX. 

35HDYPRB/ / ) 

51  FnRMAT(UX,12(F9.2,lX)/llX,5(F9.2,lX)/) 

52  F(JRMAT(11X.12(F9.2,1XH 
1  FORMAT (8F10.0) 

CALL  MASK(O) 

C  _  _ _ _ 

C  INPUT  VARIABLES 

C 

97  READ ( l >1 JHHH.TTT  »  T  ETHM. GR»  R , DT . DT 1 1  AN 
IF l HHH-1. ) 98.99.98 
99  CALL  EXIT 

_ ?8_REJtD(  1 »  1  )_PL,  T  .  1  E  TH  »D TS  » TS ,  Z_ _ _  _ 

READ ( 1 . 1  )CDC,  COB,  VOLO, MH.WBS, WDHSL ,  THE,  THED~ 

N=AN 

READ(  It  1)  (PHIOl I ),I=1,N) 

REAO( ltl)(PHI(I)fI~lyN) 

REA0(1,1)(TT(I>, 1=1,8) 

_ READ(1,1)(LLCD(I), 1=1,8) _ 
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_  01/19/68  FORTMAIN  _  _  _ 

_  READ(1,1HZZ<I)  ,1=1,24) _ _  _  _ 

_ _.REAO( 1,1)(VVG(I) ,1=1,24) _ 

_ REAO(  1 , 1 )  (  LL(  I )  ,  I  =1  ,321 _  _  _ 

. .  REAO(  1,1)  (DD(  I)  ,1  =  1,32) _  .  _  _ 

_ HEAD  (1.1)  (HWPUUI)  .1=1.321 _ _ 

II _ dZZ~ 1771  71  7 7  IN ITiALlFATlOM  _  ~  "  1.1777  177  7 771 

_ c _  .  ._  _  _  __ 

_  JJ=1  ...  ......  _ _ _  ... . . 

_  JJJ=1  _  _ _ 

_ C2=0^ _ _ 

_  XD=0. 0  _ 

_ ZD=0 .0  _  _ _ 

_ _  IBD=0.0_.  .  _ _  _ _ ..  _  _ _ 

_ THEDD=0.0 _  _  _ _ _ 

_ .  _  _  _  DO  12  1=1,  N 

_ MC  D  ( I  )=0.0 _ . _ 

_ 12  PHIDD(I)=0.’6  _  _ _ 

_ _  I  Bl-(  •.6667*WBS+PL/GR)*(».75.*V.01.Q/3,. 14159  )**  *6667  ... 

_ _ DIST  =  (.75*VCL0/3. 14159)**. 3333 _ 

_ LH=D  1ST.  _ _  _ 

_ LH1  =  LH  _  _  _ __  _  _  _ 

_ _ L  HD=0 .0 _ _ _ 

C  OUTPUT  "PROCEDURES  "7  1*7”  7" 

_ c _  .  _  .  _ _  __  _ 

_  CALL  HEAD(N)  _  _ 

. . . WRITE (3, 50.)  .  _ _ _ _ 

_ LOO_G=GR*R*R/l  Jjt+Z)  *  ( R+ZJ  ) _ 

_  _  _  IF ( Z-HHH )31 » 31, 32  _  _  _ 

.  _ _ 31  IF(T-TTT)33, 32,32  _  _ _ _ 

_ 33  I F l  TETH-TETHM)  34, 34, 32  __ _ 

_ 32  THEDEG=THE*57.  2958 _ 

THEDDE=THED*57 .2958  _ 

_ TENSJ  N)=_SQRT  (.( B-  P.Lr  Wl»Z.ppl*  L8rP  LtM  1  *ZD  PJ JJH  -  Ml  *XDP) *(H-M1»X0D)) 

_  VF(N)=-8+Ml*ZDD 


HF(N)=-H+M1*XDD 

_ IF ( N-l ) 901 ,901 ,902  _  _  _ 

902  NN=N-I  _ 

_  DO  903  1=1, NN 

_ J=N-JL _ _ 

VF(J)=-VF( J+1)-NCV(J+1)-WT(4+1)-MC(J+1)*ZCDDIJ+1) 
VF(J)=-VF(J) 

HF(J)=-HF( J+1)+NC( J+1)-MC( J+1)*XCDD( J+l) 

_  HF(J)=-HF(J)  ...  _ 

903  TENS! J)=SQRT(VF( J) *VF( J)*HF( J)*HF(j) ) 

_90l  VFW=-VF11  )-.NCV(  i)-WT(  1  )rHC(  l)*ZCPDll ) _ 

_  HFW=-HFil  )+NC(  1  )-MCU)*XCDD(  1)  _  _ 

__  TENSW=SQRT( VFW*VFW+HFW*HFW )_  _ 

WR=LCC*AN*60 • 

DO  35  1=1, N  _ _ 

__  PHIDEG(I)=PHI(I)*57.2958  _ 

_ 35  Pill  DDE  (  U  =_PH  I  0  f  I )  *57^295^. _ _ 


0  0,0  i  f  '  !  i  000 
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...  WR  ITE(  3,  36)  ..  _  ... 

36  FORMA T( /I X  ,24HRUN  ENDED  BY  CONSTRAI*  T  . / | 

WRITEC3,52) (PHIDEG(I),I=1,N) 

WRITE (3, 52) ( PHI  DDE ( I ) ,1=1 *N) 

_ WRITE(3,52)(TENS(I),I.=LiNL _  _ 

WR  ITE <  3,  51)T  »Z » X,  XD  ,  ZD,  THE.DEG  »THEODE  ,LC,*-,,  •  1.1  ENS K,  WTCT ,  H,DVB, 

IB»DYPRB  . .  .  _ 

GO  TO  97 

34  IF(JJJ-1)38, 25,38  .  _ 

38  C2=C2+1. 

_ IF(  C2-DT  1)23,24,24 _  _ _ 

24  JJ=JJ+1 

IF ( JJ-8)25,25,26 
26  CALL  HEAD(N) 

WRITE (3,50) 

JJ=1 

_25  THE  DEG =THE*57. 2958  _ _ _ _ 

_  THEDDE=THED*57.2958 

Z  Z  TENSION  IN  TETHER  ZZ  Z.Z 

ZZ  TENS(N)  =  SQRT(  ( 8-PL-M1 *ZDD J * ( B-PL-M1 *ZDD ) ♦ (H-Ml *XDD) * ( H-Ml*XDD)  ) 

_ _ VF  (  N)  =-B+  Ml *ZDD _ _  _ _ 

__  HF(N)=-H+M1*XDD 

IFIN-l ) 1 01 ,101, 102 

102  NN=N-1 

_  DO  103  1=1, NN  _  _ 

J=N- 1 

_ V FJ  J  ) =~VFJ  J ♦  1  )-NC V ( J +  1J ~WT{  J  + 1  )-MC(  J+  1  )*ZC DD  (  J+ 1 ) _ 

VF(j)=-VF(J)  _  _ 

___  HF ( J) =-HF ( J+l ) +NC ( J+l ) -MCI J+l ) *XCDD( J  +1 ) 

_  HF  ( J ) =-HF ( J ) 

103  TENS ( J )=SQRT ( VF  ( J )*VF ( J ) +HF ( J )*HF ( J ) ) 

101  VFW=- VF( l)-NCV(l)-WT(l)-MC(l)*ZCDD(l) 

_  HFW=-HF(  1)^NC(  D-MC1  1)*XCDD(1) _  _  _ 

TENSW=  SQRT ( VFW*VFW*HFW*HFW ) 

WR=LCD*AN*60. 

DO  37  1=1, N 

PHIDEG(I)=PHI(  I)*57.2958 
J37  PHIODE ( I )  =  PH  ID ( I )* 57. 2958 

_ WRITE(3,52)  (PHIOEG(I),  IM,N) _ _ _ 

WRITE (3,52) (PHI  DDE  II), 1=1, N) 

WRITE(3»52) (TENS( I ) » 1  =  1 » N) 

WRITE (3,51 )T ,Z,X,XD, ZO, THEDEG, THEDDE , LC ,WR , VG, L H, TEN SW, WTCT, H , DVB , 
1B,0YPRB 
C2=0  .0 


RUNGE-KUTTA  INTEGRATION 


23  00  74  J=  l, 4 
CALL  SUBRB 

BB(J)=FTHEDD(THED) *DT 
_  GO_J  0_  Ulx  7  2 , 7  3 17  4 ) ,  J 


.01/19/68  _  FCRTMAIM _ 

_ 71  DO  91  1  =  1,  N  _ __ _  . 

_ IR=  I _ _  _ _ 

_ CALL-SUB RC  _ _ 

_ 91 _  AAI 1 »I ) =FPHI DD (PH10)  *DT _  .... 

_ D0L81_1=1*N. _ 

_ PHI  M)=PHI  ( I  )+PHID(I)*J)T7-2. _  _ 

_  81_PHI0U.)  =  PHI0(I)+AAU, I.)Z2^ _  _ 

_ THE=THE+THED*DT/2.  ...  _ _  _ 

__  THED=THED+BBU)/2.  _  _  _  ...  . 

_  THED=O.C  _ _  _ 

_ THE=_0*0 _ 

__  _ _ T=  T+OT/2  •  ...  _  ....... 

_  DO  61  1=1, N  _  .  _  _ 

_ 6i  pH!onm=AA(i, ij/.ot _  _  .  . 

_  60  TO  74  . 

_ 72  DO  92  1=1,  N  . .  .  ... 

_ _ _ 1R=1 _ 

_ CALL  SUBRC  _  . 

_ 92_AA(2,I  )=FPHIDD(PHID)*OT _ 

_ _D0  82  1=1, N  _ _ _ 

_ PHlD(I)=PHID(I)-AAIl,n/2.4-.AA(2,J>/2. _ 

82  PHim=PHIII)+AA(  l,I|*DT/4. 

_ THE0=THE0-8Bm/2.+BB(  21/2. _ 1 _ 

_ THED=0«0  _  _ _ _ 

_  __THE=THE+BB(  1  )*DT/4# _ _ _ 

_ THE=0.0  _ _ _ _ _ 

_ DO  62  1=1,  N  _  _ _ _ _  _ 

_ 62_PHIDD(I )=AA(2,I) /CT.  _  _ 

_ GO  T0_74. _ _ 

_ 73  DO  93  1*1  VN  _ _ _ _ _ 

_ ....  IR=I  _  _ _ 

_ CALL  SUBRC  _ _ _ _ 

_ 93  AA(3,I )=FPHI00(PHI0)*DT _ _ _ 

_  _D0  83_I?C,N  _ _ _ _ 

_ PHIUJjlPHI  (U  +  DT *_{  p H.I  DlJUL/iL^rAAli ,.L1/.^±A4( 2 , JL 

_  83  PHID(  I  )  =  PHIO  ( I  )-AA (  2,1 ) /2,+AAI  3,1 )_  _ _ _  . 

THE=THE+DT*(THE0/2.-BBt  U/4. ♦B8(2)M.) 

_ THE=0. 0 

THED=THED-BB(2)/2.+BB(3) . 

_THED=0.0 

_ T=T+DT/2. _ 

DO  63  1  =  1,  N 
63  PHIDD(I)  =  AA(3,  I)/DT 
74  CONTINUE 

_  00  94  1=1, N  _  _ 

_  IR=  I  _  _ _ _ 

_ CALL  SUBRC,  _ _ 

94  AA(4, i)=FPHIDD(PHID)*OT 
DO  84  1=1, N 

PHI0(I)=PHID(I  )-AA(  3,1) 

PHI (I)=PHI ( I ) -DT  * (PHIO ( I )  +  AA(2  « I )/2 • ) 

PHI  (  I )  =PHI  (  I J+PHIDI I  )*D.T+.IAA(  l,  I  )+AA(2,I)+AA(3 
_ 84.  PH  1 0  ( L)=PH  IDC I  >  +  (_A  A.LL, .  L  L±2j.*  AA12.1. 1 1+2  ,*AA  (JL..L) 


,  I ) ) 

_+AA( 


*DT/6. 

4,I).LZ6, 


O  o 
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01/19/68  __  FORTMAIN  . .  _  . . 

THED=THEO-BB (3) 

THE=THE-DT*(THED+BB(2)/2.)  . . 

THE=TFE+THEO*DT+(BB( 1 1+8B1 2) +8B( 3) )*DT/6.  _  . 

THE=0.0  _  _  _  _ 

_ T  H  EC=T  HE  D  ♦  ( 8  B  ( 1 }  +  2  .  *B&(2 )  t  2j  ♦JBBJ-2J  +  BB141J-/  6  * _ 

THED  =  0.  C  _  _  ..  .  _ _ 

00  86  1  =  1,  N  _  .  _ 

PHIDD( I )=AA( 4, I }/DT 
MCD(I)=(PC(  I)-MC1  (in/OT 
86  MC1(  I  )=MC(  I) 

L I  NEAR  01  ST.,VEL  .  ,  AnF'aCCTjOF  SAL  LOON 

x=o.o  '  _  __ '_ '  . ~~L" 

Z=0.0  _  _ _ 

—  XD=0.  0 

_ ZD=0.0 _ _  _ _ 

xdo=o.o  __ 

_ ZDD=0 .0  .  _  __  _ 

DO  85  1  =  1, N  ______ 

X=X+LC*CGS(PHI  (I)  }  __  _ 

Z=Z+LC*SIN(PHI (I))  _  _ 

_ XD  =  XC+LCO»COS(PHI(IM-LC*SIN(PHnn>*PHIDm  __  _ _ 

XDD=XDD-2. *LC0*S IN(PHI (I) )*PHID(I )— LC*COS ( PHI ( I ) )*PHiD( I)*PHID(  I )- 
ILC*SIN(PHI(I) )*PH I D  0 ( I )  _ 

ZDD=ZDD+2.*LCD*C0S(PHI( I) ) *PHID ( I )-LC*SI N (PHI (I) ) *PHID ( I ) *PHID ( I )+ 

1 LC *C0S ( PHI ( I ) ) *PHIDC( I )  _ 

85  Z0  =  ZC*LC0*SIN(PHI(in+LC*c6S(PHiil))*PHID(  I) 

_ X=X+LH*S  IN  (T_HE ) _ 

_ Z=  Z+LH*COS  i  THE ) 

XD=XD+LH*COS(THE)*THED  _ 

_ ZD=ZO-LH*SIN(THE)*THED 

XDD=XDO-LH*SIN( THE) *THED*THED+LH*COS (THE )*THEDD 
ZDD=ZDO-LH*COS(THE)*THED*THED-LH*SIN( THE)*THEDD 

_ VB=SQRT  (XD*XD+Z0*ZD) _ _ _ _ _ 

01 ST=SQRT ( X*X+Z*Z )  _________ 

IBD=  (I8-lBn/DT 
I B 1  =  1 B 
JJ  J  =  2 

IF(TENSW+100. >97,97, 104 

104  CONTINUE  _ _ _ 

if (t-ts  > id* i i , 1 1 

11  DT=OTS 
10  GO  TO  100 
END 
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xZIZZ  ZZI  HEACINGS.roOuTPUTZ  z 

SUBROUTINE  HEA01N )  I  ~~  _ _ _ 

_ 51  FGBMAT ( 1H1» 13X»  6HPHI ( 1  )/14X» 7HPHID) I )/ 14X»7HTENS( 1 ) > 

_ 52  FORMAT (1H1,13X,6HPHI(1) »4X»6HPHII2)/14X, 7HPHID( 1 ) * 3X*  7HPHI0( 2 )/14X 

_ 1,7HTENS<  1),3X,  7HTENS12U _ _  _ 

_ 53  F0RMAT( 1H 1 , 13X,6HPHI ( 1) ,4X 16HPHI t 2) ,4X ,6HPHI { 3 )/ 14X, 7HPHID) 1) t 3X, 

_  17HPHID12 )»  3X»  7HPHID) 3J/14X* 7HTENSC 1 ) »3X» 7HTENS) 2) *3X,7HTENS( 3) ) 

__5_4_FQRMATL1HJL»  1 3X«  6HPHI  t  1)  *  4Xt  6HPHI  (2)  «4X»6HPHI  (  3  L»4X ,  6H  PH1141ZL4X  t _ 

_  17HPHID<1) ,3X,7HPHI0(2)*3Xj7HPH1D<3)*3X,7HPHID)4)/14X,7HTENS( 1),3X, 

_  27HTENS)  2  3X»  7HTENS)  3JLi.3X*7HTENS(  4 ) ).  _ _ _ 

....  55  FORMAT)  1H1,  13X,  6HPHU1)  »4X,6HPHI)2)  ,4X, 6HPHI < 3 ) ,4X ,6HPHI (4 ) ,4X , 

_  16HPHI)5)/14X,7HPHID) 1) , 3X , 7HPHID 1 2 ) f3Xf 7HPHID l 3 ) , 3X ,7HPH1D(4 ) 1 3X , 

_ 27HPHI0(5)/14X,7HTEN$U)t3Xj7HTENSI2)j3X,7HTENS(3»,3X,7HTENS(4)  »3X. 

_ 3 7H TENS) 5)) _ _ _ 

_ 56  FORMA  T)lHl,'l3X,6HPHI  J  D.«4Xa6HPHI  1 2)  »4X,  6HPHI)  3  )  ,  4X  ,  6HPHI14)  i4X* _ 

_ 16HPHIJ  5)  #4Xf  6HPHI  )6) /14X,7HPHIDli)  , 3X  , 7HPHI0)  2  ) ,  3X  ,  7HPHID)  3 ) ,  3X, 

27HPH 10  )4 )  »  3X  »7HPH  10)5  ) »  3X  J.7HPHI0)  6)  /  14X»  7HTENS  l  1 )  t  3X»  7HTENS)  2),3X, 

_  37HTENS)3)»3Xf7HT£NS)4)t3Xj  7HTENS) 5) , 3X * 7HTENS ( 6 ) ) 

_ 57  .FORMAT) 1H  1, 1 3X, 6HPHI ) 1) j4Xa6HPHI )2) t4X,6HPHI l 3 ) , 4X ,6HPHI (4)» 4X, 

_ 16HPHI )  5)  1 4X» 6HPH I  ) 6)  1 4X« 6HPHl_t 7J /  14X «7HPHID)1 )  .3X,,  7HPHIDC2U3X, _ 

_ _ 27HPHID )  3 ) «  3X«  7HPHI0(4)  »3Xi7HPJHlD) 5  )»3X*7HPHI0)6)»  3X»7HPHI0)7  )/14X, 

_ 37HTENS(l),3X17HTENS(21*3X»JHTENS)3),3X,7HTENS)4)f3X,7HTENS)5),3X, 

_ 47HTENS(6)j3X,7HTENS(7U _  _ 

_ 58_F0RMAT)IH1, 13X.6HPHI  )U14Xi6HPHl)2)f4Xf6HPHH3)f4X16HPHI(4),4X, 

_  16HPHI  l  5)  ,4Xf 6HPHI  )6)  f4Xf6HPHI<7)  ,4X,6HPHI  1 8 ) / 14X ,7HPHID) 1 )  ,  3X» 

_ 2JHPH  IDJ  2  ) j  3.X j7H P (±10.13 JjJ.X, IftPill014J  »  3 X  ,7HPH  10 )  5  )r3 X ,  7HP  HI  DJ  61,3X,._. 

_  37HPHIC  ( 7).»3X»7HPHID  )8  )/14X»7HTENS  ( 1 ) »  3X»  7HTENS  (  2 ) ,  3X,  7HTENS)  3)  *  3X, 

_ 47HT_ENS_)4  )  »3X,  7HTENS )  5 )  t  3X.r  7HTJENS)  6 ) »  3X,  7HTENS17)  »3X»7HTENS)  8)1 

_  59_F0RMATl 1H1, 13X,6HPHI< 1) > 4Xt6HPHI ( 2) ,4X,6HPHI l 3) ,4X,6HPHI 14) ,4X, 

_ 16HPHH  5)  j4X »6HPHl)6) 7)  »4X»feHPHi I  6)  *4X»6HPWI \9 ) /  14X» 

27HPHID) 1), 3Xf7HPHI0(2)f 3Xt7HPHI0(3) , 3X.7HPHI014) ,3X  ,7HPHID ) 5 ) , 3X , 
___17HPHJLQl6JLt_3X  ,.7HPH  10X7.  Ij^XjJMPljlfilfil  j3  Xt7JHPH  ID)9)/_14X ,  7H.TENS_(  1J..3X  , 
47HTENSt2) ,3X, 7HTENS13) »3X *?HTENS t 4) » 3X» THT ENSt 5) « 3X»  THTENSt 6) » 3X, 
_  57HTENS(7),3X,7HTENS)8)13X,7HTENS(9) ) 

_  _  60  FORMAT) IH1, 13X, 6HPHI ( 1) ,4X, 6HPHI ) 2) , 4X,6HPHI ) 3) ,4X,6HPHI (4) ,4X» 
LfehPHii  5)  »4X* 6HPHH6)  i4X|6HPHl i  7)  »4X,6*PHH  6)  ,4X,6HPHH9)  f4X, 

27HPHI) 10)/14X»7HPHI0( 1),3X, 7HPHID)2) ,3X,7HPHID)3)T3X,7HPHI014) ,3X, 

_ 3_7HPHIO(5)j3Xf  7HPHIDl6)TJXi7HPHIDl7)J.3X,7HPHID(8)  ,3X,7HPHID(9  )>3X  , 

48HPHI0I  iC)/14X,7'HTENSm,3X,7HTENS(2)  ,  3X ,  7HTENS)  3  )  ,3X ,  7HTENS  1 4  ) , 
53X,7HTENS(5),3X,7HTENS)6),3X,7HTENS(7)t3Xf 7HTENS) 8 ) , 3X, 7HTENSI 9) , 
63X»  8HTENS 110)) 

61  FORMAT!  1H1*  13X,6HPH1  (  i)  ,4X*6HPHU2)  ,4X,6HPM1 )  3)  ,4X,6HPHl  (4)  ,4X<» 
16HPHK5)  »4X,6HPHI  (6),4X,6HPHI  (7)  ,4X,6HPHI  )8)  ,4X»6HPHI  19)  ,4X, 

_ 27HPHI)1C)«3X»7HPHI  (_1 1)/14X,7HPHI  DIJJL*  3JC  ,  7HPHI  0  12  )  ,  3  X  ,  7  HPH 1 0 )  3  )  f  3X  , 

37HPH1D14) » 3X » THPHiOi 5) »3X» 7HPW1016 ) »3X  t7ttPWlbi7 )  » 3X  » THPH1018 ) * 3X» 
47HPHID<9> ,3X, 8HPHID110) ,2 X ,SHPHI 0 ) 1 1 ) / 14X , 7HT ENS { 1 ) » 3X , 7HTENS ( 2 ) , 
53X  »  7HTENS) 3) ,3X, 7HTENS(4) , 3X, 7HTENS (5 ) ,3X,7HTENS (6 ) , 3X , 7HTENS ( 7 ) , 

_  63X » 7HTENSI  8)  ,3X,7HTENS19)  ,3X, 8HTENSU0)  ,2X,8HTENSU1) ) 

__  62  FORMAT ( 1H1 » 1 3X»6HPHI ( 1 ) ,4X» 6HPHI ( 2 ) » 4X»  6HPHI ( 3 ) »  4X» 6HPHI (4) *  4X» 
_ 16HPHH5)  ,4X,6HPH1  (6  )  ,  4X,6HPHI  )  7 )  »  4X,  6HPHI  (  8 ) ,  AX,  6HPHI  (  9)  ,4X, 
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_ 27HPHI(10)t  3X,7HPHI(m*3X,_7HPHIU2)/14X,7HPHID(l),3X,7HPHID(2),3X, 

_  37HPHIDI3)» 3X,7HPHI0(4)J3X*JHPH10(5) »3Xt7HPHlD(6),3X,7HPHID(7),3X, 

.  __47HPHID(  8 )  1 3Xt  7HPHI01 91_i3Xt3HPHlD(  101  t2X»  8HPHI  D(  11)  ,2X  ,8HPHID(12  )/ 

_ 514X»7HTENS(  l)t  3X»7HTEN_S( 2_Lt  3Xj  7HTENS  (  3)»3X»7HTENSI4)  »3X»7HTENS  (5  )  * 

_ 63X.7HTENS<6).3X.7hTENS(7n3X«7HTENSI8).3X,7HTENSt9) ,  3X,  8HI£ALS1101t- . 


n 


72X,8HTENS(U)t2X,8HTENSIJL2il _ 

_  GO  TC  .1 1  »2  »3  »4  »5»6  t7J>8f9».10jH  *12 ).  «N 

1  WR  IT  E  l  3»  51) 

RETURN 

2  WR ITE  (3i  52 ) 


1  j 


RETURN 

_ 3  WR ITE( 3*  53)  .* 

_ 

RETURN 
WR  ITE(  3t  54 ) 
RETURN_ 

WR ITEl 3t  55 ) 
RETURN _ 


6  WR  I  TE(  3»  56  J 
RETURN 

7  WRITE! 3t  57) 
RETURN 

8  WRITE (3, 58) 
RETURN 


] 


_ 9  WRITE(3,59) 

_ RETURN 
1C  WR I TE I 3»60 ) 

_  RETURN 

11 _WR I TE ( 3 } 61 ) 
RETURN 


J 


12  WRIT  E ( 3  »  62 ) 
RETURN 
ENC 


|il 

I 
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GENERALIZED.  FORCES  ACTING  ON  BALLOON 


SUBROUTINE  SUBRB _ 

_  DOUBLE  PRECISION  TtOT  _ _ ....  _  _  ...  _ 

_ DIMENSION  AA(Af12iiBei4.bPmOE£a2)»PHIDDE(  12),MCll  12),DDD(  12) _ 

_  DIMENSION  ZC (12 ) *  SC ( 12 ) »RH0C( 12) iVC( 12  )  y VGC( 12) , XCD( 12) « ZCO( 12) 

_  REAL  MC( 12),MCD(12),LH,LC»M1»K2»MH,LCD,IB,LCLC,MCRD,MCR,LCLH, 

lLCLCO,LLCD(8),LL(32)»MMHa)jMM2I8),LCSL,NCR,NC(  12) , IBD , IB  1 ,MC 1 , LS 

_ 2.i NC  V  ( JL2L»  NQ.V R _ _  _  _  . 

_  .  COMMONTy SI »S2»S3»S4,S5»S6»S7,S8»S9»S10»SllfS12fSl3fS14»S15»S16»S17 

_ U$18,SlS,S20,S21yS22yS23.yS24t$SlySS2,S$3,SS4,SS5,PHlD(12),PHl(12), 

2THEyPFIOO( 12 ) y PHIRy  SBy DYPRB  yCDByHyB » VOLB  »RH0»G»W8»CM,CMQ,CMS»VG» 

_  3VB,THED,FTHE,FPHI,AN,AR,IJUX, Z,XDiZD,N,TT(8),RREN( 8), PHIROy VQLOy 

_  4CCD (8),VVG(24),ZZ(24) ,WDH$L»WBSyDQ(32 ) yNC ySCy VGC,VC,NCV,NCVR, 

_  5  MC , MCC»LHyLC , Ml ,  M2 , MH»LCD.  ia,LCLC»MCRD«MCR yLCL H, LC LCD  , W  WPJJU  221,  .„. 

6L  LCO, LL,MM1,MM2,TETH,S25,526,WH,WT( 12 ) yNCRyD ( 12 ) ,CDC , WTCT , TENS ( 1 2 ) 

_  7jyTENStaiDVB,PL,  ZCOC(  12  )  ,  DVC _ "  " 

_ RhOC=. 002377  _  _ 

_  CALL  CENS(Z,PR,RHG,VS) _  _ 

_  VCLB=VCL0*RHC0/RHG  _ 

_ IB=(  . 6667*hB S+PL )  *LH»LH/G _  . 

_LH  =  ( . 75* VOLB /3. 14159) **. 3333 _ 

1=1  _ 

_  512  IF ( Z-ZZ ( I) )510, 511,511 
511  1=1+1 

GO  TO  512 

_510_2SL  =,{j^Z  Zj.Zr.1.). )  /  ( Z2  ( JJ^LZ.CJLrJLU- _ _  _  - 

_  VG=VVG( I-1)+(VVG{I )— VVG (I— 1))*ZSL 

_ SB=3.1415S*(3.*V0LB/<4.*3.14159))**.6667 

XDB=XC-VG 

ZCB=ZC  _  _  ______ 

ALP=ARGD(ZDB,XDB)/57.2958 

_ CYPRB=.5*RH0*( XOB»XDB+ZOB»ZDB)  ~ _ _  _ _ _ 

H=-DYPRB*SB*LDB*COS(ALP_) . . . 

HH=WDHSL*VOLO 

hB=Hh+WBS+PL 

M1=WB/32.174+.5*RHOO*VOLO 

M2=M 1  < 

_  LS=VCL0*RH00*32.174  _  „ 

CALL  TETHER 

DVB=CYPRB*SU*CDB*SIN(ALP) 

B=L  S-HH-WB  S-DVB 

FTFE=LHMH*COS(THE)-E*SIN(THE)) 

CALL  SUMS 

_ .RETURN _ _  .... 

END 
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. .  GENERALIZED  FORCES  ACTING  ON  LINKS 

.SUBROUT  INE_SLLBRC _ ; _  _ _ 

DOUBLE  PRECISION  T,DT _ _  _ _ 

DIMENSION  AA  (4*12)  »BB(41»  PH IDEGl 12 ) » PHI  DDE  ( 12 ) , MCI ( 12) ,  DDD.M2.) 
DIMENSION  ZC  ( 12) » SC112.)  jRH0C112)  j  VC  (12)  ,VGC  ( 12  ) ,  XCD(  12) ,  2CD(  12  ) 
DIMENSION  VF  f 12) »  HF.(  121 _  _ 

REAL  MC(12) ,MCD(12).LH,LC,M1,M2,MH,LCD, IB,LCLC,MCRD,MCR,LCLH, 

1 L  CL  CD » L  L  CD.(  8 )  *  L  L.(  32J.tlLMJ.ULLi.MM2i  8 )  *LC  SL ,  NC  R ,  NC  C 1 2 )  ,  I BD  ,  LfiUHC  1 ,  LS  j 
2  *NCV ( 12 ) t  NCVR  _ 

COMMONT  »S1»S2»S3,S4»S5»S6»S7,S8»S9»S10»S11»S12»S13»S14»S15»S16,S17 
1,S18,S19,S20,S21,S22,S23,S24,SS1,SS2,SS3,SS4,SS5,PHID( 12), PHI ( 12), 
2THE,PHIDD( 12)  »PH l R , SBtDYPRB  »CDB ,H,B , VCLB, RHO, G,rfB,  CM,  CMQ, CMS,  VG, 
3VB,THED, FTHE,FPHI ,AN,ARj IR,X,2, KD, 2D,N, TT ( 6) »RREN( 0) ,PHIRD ,VOLO, 

.4CC O.ULL, V V G { 24  )^ZZ(241xMJ2HSL,ilBS,DDt-32.),NC,  SC,  VGC,VC,NCV,NCVR r~r~ 
5MC,MCD,LH,LC,M1,M2»MH,LC0»IB,LCIC ,MCRD,MCR,LCLH,LCLCD,WWPUL(32) , 
6LLCD, LL, MM1,MM2, TETH,S25»  S26, WH,HT( 12 ) ,NCR,D ( 12)  ,CDC, WTCT ,TENS (12  ) 
7,TENSW,CVB,PL,ZCDD(12),CVC,XCDD(  12) 

PHIR  =  PHI(  IR) 

PHI RD=PH ID ( I R) 

N£.R=N.C1.1R  1 _ .  .  _  _ 

NC  VR=NCV( I R )  .. _ 

MCR=MC(IR)  _ 

MCRD=MCD(IR)  _ 

CALL  SUMC  _  . 

CALL  SUMB 

.FPHl5LC*C0S(PHIR)_*(r,5*flCa*G.tBrPL-G*S25)+LC*SIN(PHIR)*(-.5*NCR-f'  .... 
1SIN(PHIR)-H)-LC«‘S26-LC*C0S(PHIR)*CQS  l PHI R ) * ,5*NCVR 

RETURN  _  _  _  _ _* _ 

END 
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SUMMATIONS. NEEDED  FOR  THETA  EQUATION 

SUBROUTINE  SUMB  _ _ 

DOUBLE  PRECISION  T,DT _ __ 

.  DIMENSION  AAI4,  12)  »B8(4).,PHIDEG(12)  » PHI  DDE!  12 ) »MC1 ( 12  )  ,DDD(  12  J 
DIMENSION  ZC 112), SCI  12) ,RHQCU2 ) ,  VC  1 12 ) ,  VGC  1 12 ) ,XCD( 1 2 ) , ZCD(  1  2 ) 
DIMENSICN  VF I 12 ) ,HF 112) _ 

REAL  MCI  12) , MCDI 12) ,LH,LC *M1, M2»MH»LCD» I B,LCLC»MCRD»MCR,LCLHi 
JLLCLCD * LLCQI 8 ) J LL L3  2 L8.) »LCSL»NCR,NC  ( 12)  ,180,1  B  L,MC  UL.S.. 

2, NCVI 12)  ,NCVR  _ _ 

COMMONT, SI ,S2»S3,S4 , S5,S6, S7, S8 , S9, S10,S11,S12,S13,SIA, S15, S16,S17 
1 , S18, S19, S20, S21 »S22  »  S23 ,S24 ,SSl , SS2  » SS3 , SS4,SS5 » PHID ( 1 2 ) , PHI ( 1 2 ) » 
2THE, PHIDD ( 12 ) ,PHI R, $B,D YPRB ,COB ,H ,8 » VOLB ,RHO,G,WB,CM,CMQ,CMS,VG, 
3VB,THED,FTHE,FPHI,AN,AR, IR, X, Z, XD, ZD, N, TTI 8 ), RREN ( 8) ,PHIRD, VOLC, 

3£C DU LlV VGJ  2AJ J,  ZZJlALtWQHil ,  OLD  132  UNC ,  S C ,  VGC » V C ,N.CV,NCV.R _+ 
5MC,M(;0,LH,LC,M1,M2,MH,LC0,  IB,LCLC,MCRD,MCR,LCLH,LCLCD,WWPUL(32)  , 
6LLC0, LL, MM1 ,MM2,T  ETH, S25, S26, WH,WT( 12),NCR,D( 12) ,CDC ,KTCT,TENS(12) 

7  ,TENSWfDVB, PL, ZCDDI 12 ),QVC,XCDD112) 

S1=0.0 

S2=0.0_  _ _ _  _ 

$3=^0 _ _  _ _ _ _ 

SA=0.0  _.  ______ 

S5=0.0_  _ 

S6=0.0  _  _  _ 

DO  1  1=1, N 

S1  =  S1+PH ID ( I ) *PH ICI I )*COS I  PHI ( I )  +THE ) 

S3=S3+PHID!H»SINIPHI  <I)+THE) _ _ _ 

S4=S4*PHI D  C I ) *COS I  PH  I II}* THE) 

S5=S5-*-SI  N  (  PHI  ( I ) +THE)  _  _  _ 

1  S6=S6+C0S ( PH  I ( I ) +THE ) 

RETURN  _ 

END  . . 
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.  SUMMATIONS  DEEDED  FOR  PHI  EQUATION 

-SUaRQim  N£.  SUK(L _ _  _ _ _ 

DOUBLE  PRECISION. T,OT _  _  ...  _  . 

DIMENSION  AA{4,12) jBB(4) ,PHIDEG( 12) ,PHIDDE ( 12 ) , MCI ( 12 ) , ODD! 12) 
DIMENSION  ZCU2),SCU2)*RHOC(  12).  VC  112)  ,VGC(  12)  ,  XCDl  12 )  ,ZCD(12  ) 
DIMENSION  VF112)*HFH2) _  _ _ 

REAL  MC( 12) »MCDI 12) *LH, LC,M1,M2,MH,LCD,IB,LCLC,MCRD,MCR#LCLH, 

1  LCLC.D » LLCD 1 8 )  «.LLX3.ZLiM(iLL6-ijilM2.(.B.)  «.LCSL«NCR  »NC  1 12),IBD,IB3*MC1,LS.. 
2.NCVI  12)  ,NCVR  _  . 

COMM0NT,SltS2»S3fS4*S5fS6jS7*S8»S9,S10tSll,S12,S13,S14,S15,S16,S17 
l.S18,S19,S20,S2lfS22,.S23iS24tSSl,SS2,SS3,SS4tSS5,PHID{121,PHI(12), 
2THE,PHIDD(12),PHIR,SBtDYPRB»CDB,HfBfVCLB,RHO,G*W9,CM,CMQtCMS,VG, 
3VBtTHED»  FTHE»  FPHI  «  AN ,  AR_j  IR  »?!jZfXD»ZD»N»TT(  8 )» RRENI  8)  ,PHIRD,VOLC  , 

.  4C.CQJL  8).,VYG(24  )tZ.L(^4iJJtfDliiLtiaStJ)i)i321JNC  ,SC,  VGC,  VC,NCV,NCVR, . . 

5MC,MCD,LH#LC,M1,M2,MH,L£D,  IBfLCLC»MCRD*MCR*LCLH»LCLCD» WWPUL (  32)  , 
6LLC0,LL,MM1,MM2,TETH,$25,S26,HH,WTC 12),NCR,D( 12) »CDC  » WTCT  »TENSl 12) 
7 ,TENSK»DVB,PL» ZCDD ( 1 2 ) »  DV C»  XCDD ( 1 2 ) 

S7=0.0 

S8=0.0  _  _ 

_ _ _ _ _  _ _  _  .  . 

510=0.0  _ 

S11=0.0 

S12=0.0  .  _ 

S  13=0.0  __  _  _ 

S14=0.0  _ 

_._$15=0.,Q _  .  _  _ 

S16=0.0 

S 17=0 .0  _  __  _ 

SI  8=0 .0 
S 19=0 .0 
S20=0.0 

_ S  2.1=0  *0 _ _ .......  .  _ _ _ _ 

S22=0.0 
S23=0.0 
S25=0 .0 
S26=0.0 
SS1=0.0 

SS2=0.0 _ _ _  _ _ 

SS3=0 .0 
SS4=Q.O 
SS5=0 .0 
DO  1  1=1,  N 

S7=S7+PHI0( I )*PHIG(I)*SIN(PHIl I)+THE) 

S8=S6+PHIDm»C0S  (PHI  (I)-PHIR)  _ 

S22=S22*PHID l I )*PHID( I ) ♦SIN! PHI ( I ) -PHI R) 

1  CONTINUE 

IFU-IR+1  )9,9, 11 
9  J=I R-  1 
DO  2  1=1, J 

_ S9=S9+PHI0(  I  )»COS  (PHU  D-PHIR) _ 


I 
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_ S11=S11+SIN(PHI  ( I  J-P.HI8J _ _ _  _  . 

S12=S12+PHIDU)*PHI0jm*SIN(PHlU)-PHIR) 

„  2  CONTINUE  . 

11  CONTINUE  _  _ 

_ IFJJR  +  l-N  )  8»  8_t  6 _  _  _ _ _ 

_  8  J=I R+ 1  _ 

00  3  I=J,  N  _  _ _ _  _ 

_ I N=I- 1  _  _  _  . 

DG  7  11=1, IN 

SSl=SSl+MCOU)*PHlD(in*COS(PHI  (In-PHIR) 

_ SS_2=S S 2+MCJJ  )*PHIDUI)»CQS(PHI(I I  )-PH I R )  _  . 

SS  3=SS3+MCD{ J)*SIN{PHI(II )-PHIR ) 

_ 7  SS4=SS4+MC(J)*PHID(II  )*PHIDUI)*SIN(PHI(I  IJ-PHIR) 

__  3  CONTINUE  _ _ 

J  =  IR  + 1  .  _ _  _ _ 

00  4  IN=  J  , N 

_ S 17=$  1.7+MCO  (  IN )  *PHJ  DI1N)*C0SI~PHH  IN  1_-_PH  IR)  .  _  _ _ 

S18=S 18+MCl IN)*PHID( IN)*COS(PHI ( IN J-PHIR ) 
S20=S20+MCD(IN)*SIN(PHI1 INJ-PHIR) 

■_  S21r$21+MC(  IN)*PHID(  IN)*PHID{  INj*SIN{PHI(IN)-PHIR) 

S25=S25+  MC ( I N)  _  .  _  .. 

S26=$26+  (  NC  ( I N)  *S  I NI  PHI  U  N) )  *S  I N  { PHI  {  IN  )  )  +NCV  l  IN  )*COS  (  PHI!  IN)  )* 

_ ICOStPHU  IN ) ) ) *C0 S ( PHI R-PH I C  IN ) ) _ _  .....  _  _ 

_ 4  CONTINUE 

fc  CONTINUE  _  _ 

RETURN  _  .  .  _ 

END 


|noo  i  I  .  ooo 
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_  FQRCES-.ACTI NG  ON  TETHER 

_ SUB  BOUT  l  NE_TETHE8 _ _ _ _ 1 _  _ 

_  DOUBLE  PRECISION  T,OT _ _  .  _ _ 

DIMENSION  AA(4,121,BB(A1  iPHIOEGU2>  ,PHIDDE(12  > , MCI  1 12  J  ,DDDU2) 
DIMENSION  ZC( 12) ,  SCI  12) iRHQCI 12 )> VC (12) , VGCI12) , XCD< 12) , ZCD112 ) 
DIMENSION  VFU2)  ,HF(12)  _  . .  . 

REAL  MCI  12) »MCD 112) ,LH,LC,M1, M2,MH,LCD,18,LCLC,MCRD,MCR,LCLH, 

_ 1 LC1£D , L LC DID .) J.LL.L 3 2JaJ4MllBljiiM2J8_l  ,JL C_SL» NC R * NC 1 12) » I BDj LB 

_ 2*NCVl  12)  f  NCVR  _ _ 

COMMONT, S1,S2,S3, S4»S5,S6,S7, $8, S9, SIO, S 1 1, S 12,S 13, S14, S 15, S16, S 17 
1,S18,S19,S20,S21, S22,S23,S24,SS1,SS2,SS3,SS4,SS5,PHID(12),PHI(12), 
_  2THE,  PHIDDI 12 ) , PHIR , SB , DYPRB ,CDB * H ,B ,VOL8,RHO,G,W6,CM,CMU,CMS, VG, 

.  3VB»THE0»FTHE,FPHI,AN,  ARtURjX*.  Zt.XDi  ZD»N*TT{  8 ) » RREN I  8),PHIRD»V0LC» 

_ 4CC  D.L8  Lt  VVJG1 2  4J  ,Z  U  24JL*HQliSi.  j  WBiiDD  L32  U  NC , S  C ,  V  GC ,  V  C ,  N  CV ,  N  CVR , _ 

5MC,MCD,LH,LC,M1,M2,MH,LCD, I  3, LCLC ,MCR D, MCR, LCLH, LCLCD , WWPUL  (  32) , 
6LLCD,LL,MM1,MM2,TETH,S25,S26,  WH,WTI  12  )  ,NCR ,D( 12)  ,CDC  ,  WTC  T ,  TENS  (  1 2 ) 
7,TENSW»DVB,PL,  ZCDD ( 12  )  ,  DVC,  XCDDI 12 ) 

LINK  LENGTH, RATE  OF  CHANGE, MASS , AND  REF.  AREA. 


_  1=1  _ _ 

512  I F ( T —  TT ( I ) ) 5 10 , 5 1 1 , 5 11 

511  1=1  +  1  _  _ _ 

GO  TO  512 

510  111  =  1-1  ..  .  _ 

_ TETH=Q.,0 _ _ _ 

IF! 111-2)515,514,514 

514  DO  513  11=2, III  .  _ 

513  TETH=TETH+LLCD( I  I ) *( TTI 1 1 )-TT ( I  1-1) ) 

515  Tfc  TH=TE  TH+LLCD  (  I  )  4 1  T-TT(.I  -1 )  ). 

LC=TETH/ AN  _ 

_ L  C.D=_L  L  CD  (  I.)/  AN _ _ _ 

D2=DD(  1) 

WPUL=WWPUL (  1  ) 

WTCT=0. C 
DO  599  J=1 , N 
A  J= J 

_ 1  =  1  _  _ _ 

612  IFIAJ*LC-LL( I ) ) 6 10, 6  1 1, 61 1 
611  1=1+1 

GO  TO  612 

61 C  LCSL=(LC*AJ-LL(I-1 ) )/ (LL( I ) — L L C 1-1 ) ) 

0  1=D2 

WPUL_1=WPUL_ _ _ _ 

WPUL= WWPUL  (  I  —  1 )  +  (  WWPULII  )' -WWPULlI-1  )  )  *  L  C  S  L 
D2=DD(  1-1  )  +  (DD(  I  )  —  DD C  1-1)  )*LCSL 
DDD (J)=(Dl+D2)/2. 

WT (J  ?  =  ( WPUL  +  WPUL1 ) *LC /2 . 

599  WTCT=W TCT+WT ( J > 

DO  598  1=1, N 
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__  ZC(I)=0.0 

_  xcom=o.o 

_.  ZC0(II=0.0 


_  ZCOOt  I  )=0.0  _  _ 

_ xcpoj  il=o.o _ 

_  LLL=N+1- I  _ 

_ D( LLl ) =DDD( I )  _  _ 

598  MC ( LLL) =WT( I ) /32. 1 74 
DO  100  1*1,  N 
WT ( I)=MC(I)*32.174 

_ SC(_U=LC*D(I  ) _ 

C..„ .  _  LINEAR  DYNAMICS  OF  LINKS 

IF (1-1 ) 103,103* 101 

101  I  1=  1-1  _  _ _  _ _  _ _ 


_ DQ_1Q.2_J=L.I  I _ _ 

_ XCDU)=XCD(  I  )+LCD*COS(  PHUJ.)  LtJ.C*$IN(  PH  I  ( J  )  )*PHID(  J) 

_ ZC  D  ( I )  =  Z  CD  ( I )  +L  C  D*  S  I  Mi  PH.IJ.J. >  ).+L C*  CO  S  (  P  H I  (  J  )  )  * P H I D  (  J ) 

__  ZCDD( I I=ZCDD ( I )+LCD*COS (PHI (J ) )*PHID( J) *2. +LC*COS( PHI ( J) ) *PHI DD(  J) 

_ 1— LC*SI N( PHI ( J) ) *PHID( J)*PHID( J) 

XCDD( I )=XCDC ( I ) -2  «*LCD*S IN ( PHI ( J ) )*PHID(J  )-LC*COS( P HI ( J ) ) *PHID( J ) * 
_ 1PHI01 J )-LC*SIN(PHI ( J) ) *PH1 DDt  J) _  .  _  . 

102  ZC(  I)=ZC(  I  )+LC*SIN(PHI(4n  _ 

103  ZC(I  )  =  ZC(I  )  +  .5*LC*SIN(PHUi)  ) 

XCOI.I  )=XCD(  I  )  +  .5*LCD*C0S(PHlU)  )r  -  5*LC*SI  N(  PH  I  (  I  ))*PHID(I  ) 

ZCDU  )=ZCD(  I  )  +  .5*LCD*SlN(PHl(  I )  )  + . 5*LC*C0 S( PH  I(  I  )  )*PHID(I) 

ZCDO(  I  )=  ZCDD  l  I  )+LCD*COS(PHI(  I  ) )*PHlD(I)+.  5*LC*C0S( PHI ( I ) ) *PHI  DD(  I  ) 


XCDD  (  I  )=XCDO  ( I  )-LCO*S  IN  (PHI  (I )  )*PHID.(  I )- .  5*LC*C0  S  (  PHI  (  I)  )*PHID(I)* 
IPHID( 1)~.5*LC*SIN(PHI ( I ) > *PHI DD ( I ) 

C  .  .. 

C  _  _  .  DRAG  FORCEi._0N. .LINKS 

.CAl  L_OJJiS,(  ZC  (  n,PR,RHOC(  I),VSC) _ _ _ _  _ _ _ 

100  VC(I)  =SCRT(XCD(I)  *XCD(I)*ZCDm*ZCD(I  )  ) 

DO  713  J=  1  »N 
1  =  1 

712  IF (ZC ( J)-ZZ( I ) ) 710,711,711 
711  1*1+1 

_ _  GO  T0_  712  _ 

710  ZCSL* ( ZC ( j )-ZZ ( I- 1) )/(ZZ( I )-ZZ( 1-1)1 
VGC ( J ) =VVG  f I-1)  +  (VVG(I )-VVG( 1-1) )*ZCSL 
NCV (J )=CDC*.5*RH0C( J )*SC( J)*ZCD( J  )*ZCD( J) 

713  NC ( J )=CDC* • 5*RH0C ( J ) *SC ( J ) *( ( VGC (J)-XCD(J))*( VGC ( J )-XCD( J ) ) ) 

DO  20C  J= 1 , N  _ 

_ I F I ZCC ( J ) ) 202 , 202 , 20  3 _ 

202  NCV ( J )=-NCV ( J  j  _ _ 

203  CONTINUE 

DVC=0 .0  _ _ 

DO  204  1=1, N  _  _ 

204  DVC=DVC+NCV( I) 

IF  <VGC(J)-XCD(J))201j  200J.200 _ 


l 
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I  .  _  ....  2C1  NCI.J)=-NCI  J)  _  . 

f.  .....  200  CONTINUE 

,  LCLC=LC*LC  .......  _ 

I  LCLH=LC*LH  ........  _ 

I _ L  CtC  Or_LC*LCD _ 

L  ...  .._  RETURN 

L  END 


